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ABSTRACT 


The  objective  of  this  research  is  the  performance  of  quantitative 
infrared  spectroscopic  measurements  of  boron  fluorides  and  oxides  in  the 
gas  phase.  The  measurement  program  is  designed  to  determine  infrared 
absorption  coefficients  and  absorption  line  spacings  in  order  to  facilitate 
the  construction  of  infrared  radiation  band  models. 

In  order  to  carry  out  these  measurements  an  apparatus  which  permits 
the  performance  of  quantitative  spectroscopic  measurements  of  transient 
species  over  the  temperature  range  of  300  to  1500  K  was  constructed.  The 
capability  of  this  apparatus  to  perform  these  measurements  over  its  designed 
range  of  operating  parameters  was  demonstrated. 

The  infrared  absorption  coefficients  of  the  nu-3  band  of  BF^  were 
measured  over  the  range  of  300  to  1500  K.  These  measurements  were  combined 
with  flame  measurements  at  2400  K  to  develop  and  verify  a  radiation  band 
model  for  BF^.  The  formalism  developed  for  this  band  model  is  more  comprehen¬ 
sive  than  previous  treatments  and  is  generally  applicable  to  perpendicular 
bands  of  symmetric  top  molecules. 

The  infrared  band  strength  of  BF  was  measured.  The  band  strength  was 
determined  by  measuring  the  absorption  cross  section  of  a  single  vibration- 
rotation  line  using  a  tunable  diode  laser  as  the  light  source.  The  BF  number 
density  was  measured  in-situ  using  ultraviolet  absorption  spectroscopy. 


1.  INTRODUCTION 


The  objective  of  this  research  program  is  the  performance  of  quantitative 
infrared  spectroscopic  measurements  on  a  number  of  boron  oxides,  fluorides, 
and  oxy-fluorides.  The  measurement  program  is  designed  to  produce  results 
which  are  suitable  for  use  in  infrared  radiation  band  models,  and,  therefore, 
the  infrared  absorption  measurements  may  be  performed  at  high  spectral  resolu¬ 
tion,  as  well  as  at  low  spectral  resolution.  This  permits  the  determination 
of  the  two  band  model  parameters,  the  average  absorption  coefficient  and  the 
fine  structure  parameter.  The  measurements  are  performed  within  the  temperature 
range  of  300  to  1500 K.  The  candidate  molecules  for  the  study  are  BF^,  BF,  BO2, 
OBF,  and  HBO2.  With  the  exception  of  BF^  which  is  a  commercially  available  gas, 
these  are  transient  species  which  must  be  produced  in  high  temperature  reactions. 

This  research  has  produced  several  significant  accomplishments.  These  are: 

(1)  The  development  of  a  versatile  apparatus  for  the 
study  of  the  spectroscopy  and  kinetics  of  transient 
species  between  300  and  1500 K, 

(2)  The  measurement  of  the  infrared  spectral  absorption 
coefficients  of  the  BF^  band  between  300  and 
1500 K  and  the  development  of  a  radiation  band 
model  for  this  band,  and 

(3)  The  measurement  of  the  fundamental  infrared  band 
strength  of  BF. 

The  apparatus  that  was  developed  under  this  contract  is  basicly  a  high 
temperature  flow  reactor  which  is  suitable  for  the  performance  of  quantitative 
spectral  measurements  of  transient  species  over  the  temperature  range  of  300  to 
1500 K.  In  addition  to  quantitative  infrared  spectroscopy  at  low  or  high  spectral 
resolution,  the  apparatus  is  capable  of  studying  transient  species  by  ultraviolet 
absorption  spectroscopy,  laser- induced  fluorescence,  microwave  discharge  lamp 
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induced  fluorescence,  and  molecular  beam  mass  spectrometry.  Finally,  it 
should  be  noted  that  the  apparatus  can  also  be  used  to  study  the  reaction 
kinetics  of  transient  species.  This  apparatus  is  described  in  Section  2 
which  is  a  preprint  of  an  article  that  will  be  submitted  for  publication 
in  the  Review  of  Scientific  Instruments. 

The  spectral  absorption  coefficients  of  the  BF^  band  were  measured 
over  the  temperature  range  of  300  to  1500  K.  An  advanced  infrared  radiation 
band  model,  which  is  generally  applicable  to  perpendicular  bands  of  symmetric 
top  molecules,  was  formulated  for  BF^.  The  spectral  absorption  coefficients 
measured  in  this  study,  along  with  those  measured  in  a  flame  at  2400 K,  were 
used  to  refine  and  verify  the  validity  of  the  band  model.  These  results  are 
presented  in  Section  3  which  is  a  preprint  of  an  article  that  will  be  submitted 
for  publication  in  the  Journal  of  Chemical  Physics.  A  preliminary  version  of 
this  work  appeared  previously.  (M.E.  Gersh  and  C.E.  Kolb,  "Infrared  Spectral 
Measurements  of  Species  in  a  High  Temperature  Flow  Tube,"  in  Proc.  lQth  Materials 
Research  Symp.  on  Characterization  of  High  Temperature  Vapors  and  Gases,  J.W.  Hasti{„. 
ed.,  NBS  Special  Publication  561  (1979),  p.  679.)  r 

The  fundamental  infrared  band  strength  of  BF  was  measured  using  a  tunable  , , 

diode  laser  as  the  radiation  source.  The  BF  was  created  in  a  2450  MHz  micro- 
wave  discharge  of  a  mixture  of  BF^  and  He  external  to  the  flow  reactor.  The 
band  strength  was  determined  from  a  measurement  of  absorption  cross  section  of 
the  P(ll)  line  of  BF  at  300 K.  The  BF  number  density  was  measured  in-situ  by 
ultraviolet  absorption  spectroscopy  in  the  BF  A1^  x  *£  band  head  at  196  nm. 

These  results  are  presented  in  Section  4  which  is  a  preprint  of  an  article  that 
will  be  submitted  for  publication  in  the  Journal  of  Chemical  Physics.  -  * 
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2.  A  VERSATILE  HIGH  TEMPERATURE  FLOW  REACTOR  FOR 
KINETIC  AND  SPECTROSCOPIC  STUDIES 


by 


Michael  E.  Gersh,  Joel  A.  Silver,  Mark  S.  Zahniser, 
Charles  E.  Kolb,  Robert  G.  Brown,  Conrad  M.  Gozewski, 
Spiros  Kallelis,  and  Joda  C.  Wormhoudt 


ABSTRACT 


A  high  temperature  flow  reactor  has  been  constructed  to  perforin  kinetic 
and  spectroscopic  studies  on  species  at  temperatures  between  300  and  1500K. 
The  species  in  the  flow  reactor  may  be  detected  or  studied  using:  laser 
induced  fluorescence,  resonance  fluorescence,  infrared  absorption,  ultraviolet 
absorption,  chemiluminescence,  or  molecular  beam  mass  spectrometry.  The  fea¬ 
tures  of  the  apparatus  are  described  in  detail,  and  the  initial  spectroscopic 
and  kinetic  studies  are  discussed. 
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INTRODUCTION 


Flow  Cubes  have  been  used  routinely  Co  perform  kinetic  measurements  at 

temperatures  between  200  and  600  K,  and  flame  and  shock  tube  techniques  are 

commonly  employed  for  temperatures  above  1000  K.  However,  relatively  few 

data  are  available  in  the  500  to  1500  K  regime,  which  bridges  the  gap  between 
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low  and  high  temperature  measurements. 

On  the  other  hand,  flow  tubes  have  not  generally  been  employed  in  the 
performance  of  quantitative  infrared  measurements.  These  have  most  often 
been  performed  in  static  cells,  shock  tubes,  and  burners. Static  cells 
are  not  convenient  for  the  study  of  transient  species  and  are  limited  to 
approximately  1200  K.  Shock  tubes  and  burners  are  not  generally  useful 
below  1000  K,  and  they  usually  require  equilibrium  calculations  of  species 
number  densities. 

This  paper  describes  a  high  temperature  flow  reactor  which  was  developed 
to  address  both  of  these  problems.  It  is  designed  to  perform  kinetic  and 
spectroscopic  studies  over  the  temperature  range  of  300  to  1500  K.  The  kin¬ 
etic  measurements  are  facilitated  by  the  capability  to  measure  reactant  and/or 
product  number  densities  using  a  wide  range  of  diagnostic  techniques:  laser 
induced  fluorescence,  discharge  lamp  resonance  fluorescence,  infrared  absorption, 
ultraviolet  absorption,  chemiluminescence,  and  molecular  beam  mass  spectro¬ 
metry.  In  addition,  the  apparatus  was  designed  to  perform  quantitative  infrared 
absorption  measurements,  with  the  application  of  one  or  more  of  the  other  diag¬ 
nostic  capabilities  to  independently  measure  the  number  density  of  the  absorb- 


ing  species.  As  a  corollary,  the  multiplicity  of  diagnostic  capabilities 
also  permits  the  performance  of  visible/ultraviolet  spectroscopic  measurements 
of  transient  species. 


I.  DESCRIPTION  OF  THE  APPARATUS 
A.  General 

The  high  temperature  flow  reactor  has  been  designed  with  the  capability 

to  perform  both  kinetic  and  spectroscopic  studies  of  species  at  temperatures 

between  300  K  and  1500  K.  As  a  result  the  apparatus  design  is  considerably 

more  complicated  than  would  be  required  if  only  kinetic  studies  were  to  be 
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performed.  ’  *  This  is  due  to  the  fact  that  only  relative  species  concentra¬ 
tions  need  be  measured  to  determine  (pseudo-first  order)  rate  constants,^  but 
absolute  species  concentrations  and  well-defined  path  lengths  are  required 
for  the  performance  of  quantitative  spectroscopic  measurements.  In  addition, 
the  apparatus  has  been  designed  to  accomodate  a  great  variety  of  measurement 
techniques,  thereby  providing  the  flexibility  to  study  a  wide  range  of  species. 

A  block  diagram  of  the  major  apparatus  components  is  shown  in  Fig.  1.  A 
gas  handling  system  external  to  the  flow  reactor  is  used  to  prepare  the  several 
gas  flows  that  are  introduced  into  the  reactor.  One  (or  more)  of  these  flows 
may  be  passed  through  microwave  discharges  to  produce  radical  species.  In  the 
flow  reactor,  gases  may  be  mixed  to  produce  the  reactions  whose  rates  are  to 
be  measured,  or  reactions  may  be  carried  out  in  order  to  produce  transient 
species  for  spectroscopic  studies.  The  flowing  gases  then  pass  into  an  analysis 
region  where  a  variety  of  diagnostic  techniques  may  be  employed  to  characterize 
the  species  present.  The  following  measurements  may  be  performed:  infrared 
absorption  using  either  resonance  or  arc  lamps,  ultraviolet/visible  discharge 
lamp  induced  resonance  fluorescence,  ultraviolet/visible  laser  induced  fluor¬ 
escence,  and  molecular  beam  mass  spectroscopy.  This  versatility  permits  the 
study  of  almost  any  species. 
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B.  Vacuum  Jacket  and  Flow  Reactor 
B.l.  Overall  description 

The  detailed  design  of  the  variable  temperature  flow  reactor  is  shown 
in  Fig.  2.  The  outer  vacuum  jacket  is  constructed  of  steel  and  is  water 
cooled,  as  are  the  aluminum  end  plate  assemblies.  The  water  cooling  protects 
temperature-sensitive  components  and  insures  dimensional  stability  for  optical 
components  (such  as  the  White  cell  mirrors)  mounted  on  the  vacuum  jacket  flanges 
The  upstream  end  plate  holds  the  alumina  flow  tube,  as  well  as  power,  gas,  and 
thermocouple  feedthroughs.  The  downstream  flange  is  used  for  experiments  for 
which  mass  spectrometer  sampling  is  not  required.  When  the  mass  spectrometer 
is  employed,  the  flow  reactor  is  mated  directly  to  the  outer  vacuum  wall  of 
the  mass  spectrometer.  In  the  analysis  region,  eight  10  cm  inside  diameter 
flanges  are  spaced  equally  around  the  circumference  of  the  vacuum  jacket.  Four 
of  these  ports  are  connected  by  alumina  tubes  to  the  main  flow  tube  and  provide 
for  optical  analysis  of  the  hot  flowing  gas.  One  of  the  remaining  four  ports 
is  presently  connected  to  a  Kinney  KMBD  1602  mechanical  pump  and  Roots  blower 
with  an  effective  pumping  speed  of  approximately  1000  cfm  (450  £/sec) . 

The  main  elements  of  the  thermal  design  include  the  heating  elements 
and  the  surrounding  layer  of  refractory  insulation.  The  insulation  is  composed 
of  a  self-supporting  "spool"  of  rigid  alumina  fiber  insulation  (Zircar  Products, 
Inc.),  upon  which  is  wound  several  layers  of  Zircar  alumina  fiber  blanket  insula 
tlon,  followed  by  additional  layers  of  alumina-silica  blanket  insulation 
(Carborundum  Corp.).  The  insulation  assembly  has  an  inside  diameter  of  15  cm 
and  Is  designed  to  withstand  continuous  operating  temperatures  of  1800  K.  The 


commercially  available  Kanchal  A-l  heating  elements  (Thermcraft,  Inc.,  type  RH) 
have  a  design  temperature  rating  of  1500  K  and  are  operated  at  120  VAC.  The 
heating  elements  are  arranged  in  a  three  zone  configuration  to  minimize  tempera¬ 
ture  gradients  due  to  end  effects.  The  heated  length  is  78  cm,  with  a  total 
power  capability  of  approximately  8  kW.  The  heaters  are  supplied  with  power 
by  solid-state  proportional  temperature  controllers  and  power  supply  units 
(Love  Controls,  Corp.).  The  fully  assembled  flow  reactor  achieves  a  tempera¬ 
ture  of  1500  K  with  a  power  consumption  of  approximately  3  kW.  The  thermal 
stresses  developed  in  the  alumina  flow  tube  components  are  small  enough  that 
the  apparatus  is  routinely  heated  from  room  temperature  to  1500  K  in  approximately 
4-6  hrs,  and  no  damage  to  the  alumina  components  has  occurred  during  the  ll-$ 
years  that  the  heated  flow  reactor  has  been  in  use. 

The  design  of  the  gas  flow  system  within  the  variable  temperature  flow 
reactor  is  relatively  complicated  due  to  the  constraints  required  to  provide 
the  capability  of  performing  quantitative  spectroscopic  measurements  in  the 
analysis  region.  The  main  flow  in  the  7.3  cm  inside  diameter  99.8%  alumina 
flow  tube  is  usually  composed  of  He  or  Ar  carrier  gas.  This  gas  is  heated  as 
it  passes  through  two  heat  exchangers  which  are  2.5  cm  thick  mullite  "honeycomb" 
hole  arrays  with  1.6  mm  round  holes  and  a  67%  frontal  open  area  (General 
Refractories  Co.).  These  heat  exchangers  are  required  to  heat  the  carrier  gas 
at  high  flow  rates  (large  Reynolds  numbers) ,  and  they  also  serve  to  laminarize 
the  carrier  gas  flow.  Downstream  of  the  heat  exchangers,  two  concentric  99.8% 
alumina  injectors  permit  the  introduction  of  reactive  species  into  the  flow. 
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The  gas  in  the  Inner  injector  (9.5  mm  i.d.,  12.7  mm  o.d.)  may  be  passed  through 
a  commercial  2450  MHz  microwave  discharge  cavity  (Ophthos  Instrument  Co.)  before 
being  introduced  approximately  1-2  cm  upstream  of  the  end  of  the  outer  injector 
(14.3  mm  i.d.,  19.1  mm  o.d.).  The  microwave  discharge  is  used  to  generate  atoms 
and  free  radicals  which  can  be  injected  directly  into  the  main  carrier  gas  flow 
or  can  be  mixed  with  a  reactant  gas  in  the  outer  coaxial  tube  before  being  injected 
into  the  main  flow.  This  arrangement  permits  a  great  deal  of  flexibility  in  the 
preparation  of  reactant  species  for  study.  The  heated  gas  then  passes  through 
an  isothermal  region  of  the  flow  tube  with  constant  cross  section.  In  this 
region  a  second  reactant  species  may  be  introduced  into  the  flow  via  a  movable 
injector  for  the  measurement  of  rate  constants.1  The  present  movable  injector 
is  constructed  of  6.3  mm  o.d.  quartz  tubing  (with  a  temperature  limit  of  approxi¬ 
mately  1250  K).  The  gas  is  injected  into  the  main  flow  from  a  5  cm  o.d.  loop  of 
tubing  having  16  0.4-0. 5  mm  diameter  holes  to  promote  rapid  and  uniform  mixing. 

The  main  flow  with  its  mixture  of  reactive  species  then  passes  through 
a  subsonic  nozzle  (area  ratio  3.3:1,  4.0  cm  i.d.)  to  enter  the  analysis  region, 
as  shown  in  Fig.  3.  A  second,  external  flow  (shield  flow)  of  pure  carrier  gas 
is  heated  to  the  temperature  of  the  main  flow  as  it  passes  between  the  heating 

elements  and  the  outside  of  the  flow  tube.  Then,  after  the  shield  flow  passes 

through  a  subsonic  nozzle  which  is  concentric  with  the  main  flow  nozzle,  the 

flows  merge.  The  area  ratio  of  the  main  flow  nozzle  to  the  shield  flow  nozzle 

is  1:1,  and  the  corresponding  flow  rates  should  ideally  also  be  equal  so  that 
the  velocities  of  the  two  flows  will  be  the  same  when  they  meed  downstream  of 
the  nozzles  (in  order  to  minimize  mixing  due  to  shear  forces).  This  results 


in  a  fluid  dynamic  confinement  of  the  reactive  species  in  the  main  flow  as 
they  pass  through  the  analysis  region.  In  addition,  small  flows  of  carrier 
gas  are  added  to  each  of  four  optics  tubes  as  a  purge  gas  in  order  to  keep 
the  optics  clean.  The  flow  of  purge  gas  into  the  optics  tubes  is  generally 
set  at  approximately  10  percent  of  the  main  flow  (equally  divided  between  the 
four  optics  tubes),  although  this  ratio  is  not  critical.  The  six  flows  of 
gas  meet  in  the  analysis  region  volume  which  is  defined  by  a  six-arm  cross 
piece  that  is  made  from  machinable  alumina  ceramic  (Cotronics  Corp.  type  960). 
The  cross  piece  supports  the  four  optics  tubes  and  the  flow  tube.  (The  front 
of  the  flow  tube  rests  on  two  sapphire  balls,  thereby  permitting  the  flow 
tube  to  expand  without  generating  any  stress  on  the  cross  piece.)  However, 
the  main  purpose  of  the  cross  piece  is  to  provide  an  environment  in  which  the 
main  flow  from  the  flow  tube  will  be  well  behaved  as  it  traverses  the  optical 
analysis  region,  impinges  on  the  mass  spectrometer  sampling  orifice,  and  is 
exhausted  by  the  mechanical  pump.  Thus,  the  design  of  this  element  of  the  flow 
system  is  particularly  important  to  the  success  of  the  experimental  program, 
and  separate  experimental  studies,  described  below,  were  undertaken  in  order 
to  insure  that  the  design  would  adequately  control  the  gas  flow  in  the  analysis 
region. 

B.2.  Flow  simulation  experiment 

A  crucial  element  of  the  apparatus  design  is  the  configuration  of  the 
flow  system  in  the  spectroscopic  analysis  region.  In  order  for  one  to 
be  able  to  perform  quantitative  spectroscopic  measurements  on  the  species 


2-8 


in  the  main  flow,  three  criteria  must  be  satisfied:  (1)  geometry  of 
the  gas  flow  must  be  well  defined  as  it  traverses  the  analysis  region,  (2) 
the  flow  in  this  region  must  be  stable,  and  (3)  the  temperature  of  the  gas 
flow  must  be  well  defined  (preferably  isothermal).  These  criteria  may  be 
satisfied  by  maintaining  a  stable  cylindrical  flow,  without  recirculation 
of  absorbing  species  in  the  optics  tribes.  For  high  temperature  gas  flows 
and/or  multi-pass  absorption  cells,  this  necessitates  the  use  of  a  wall-less 
confinement  of  the  absorbing  gas  species  as  they  pass  through  the  analysis 
region.  This  can  only  be  accomplished  by  fluid  mechanical  confinement  of 
the  flowing  gas. 

Since  this  portion  of  the  apparatus  represented  a  novel  design,  it  was 
felt  that  its  design  should  be  tested  in  a  flow  simulation  experiment  before 
the  fabrication  of  the  six-arm  cross  piece  out  of  machinable  alumina  ceramic. 
This  obviated  the  need  for  a  major  system  redesign  if  the  initial  design  proved 
to  be  inadequate  and  would  establish  the  approximate  range  of  operating  flow 
conditions  over  which  the  fluid  mechanical  confinement  of  the  absorbing  species 
would  be  effective. 

t 

The  flow  simulation  experiment  was  set  up  in  the  following  manner:  a 

half  scale  model  of  the  flow  system  was  fabricated  out  of  plastic  (polymethvl 

methacrylate)  with  flowmeters  monitoring  the  main  flow,  the  shield  flow,  and 

the  four  purge  flows.  The  gas  flow  was  simulated  by  using  a  water  flow  with 

the  same  Reynold's  number,  which  is  a  standard  technique  in  fluid  mechanics  and 

9 

is  valid  for  subsonic,  incompressible  flows.  The  flow  pattern  is  determined 
by  injecting  Fluorescein  dye  into  any  one  of  the  water  flows  and  observing 
the  bright  fluorescence  induced  by  illumination  with  a  tungsten/halogen  lamp. 
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The  water  flow  velocity  required  for  the  simulation  is  obtained  from  the 
flow  velocity  of  (for  example)  Ar  that  is  to  be  simulated  by  equating  the 
Reynolds  numbers  of  the  flows: 

R  e  .  .  HXil  (1) 

n  n' 

where  p  is  the  density,  u  is  the  flow  velocity,  d  is  the  tube  diameter, 

n  is  the  viscosity,  and  the  prime  refers  to  the  conditions  for  water.  Then, 

we  let  d"  =  d/2  and  p  -  pq  P,  where  pQ  -  2.35  x  10~6  g/cm3/Xorr  and  P  is  the 

Ar  pressure  in  Torr.  Also,  n'  -  1.0  g/cm3,  n*  2.2  x  10"4  poise,  and 
_2 

n'  =  1.0  x  10  poise  (at  293  K) .  This  leads  to  the  expression  for  the  water 
flow  velocity  u'  required  for  the  simulation: 

u'  *  2.1  x  10  Su  (2) 

To  simulate  Ar  flows  at  other  temperatures,  one  would  have  to  substitute  the 
appropriate  values  of  po  and  n. 

The  flow  simulation  experiment  was  actually  performed  for  a  slightlv 
different  apparatus  configuration  from  that  shown  in  Fig.  3.  In  particular, 
the  nozzle  for  the  main  flow  had  an  area  ratio  of  2.1  (ratio  of  main  flow 
to  shield  flow  *  2:1), 10  as  opposed  to  the  area  ratio  of  3.3  currently  usee 
in  the  high  temperature  flow  reactor.  The  flow  reactor  was  initially  constructe 
with  the  smaller  area  ratio  nozzle,  but  this  design  proved  to  be  marginal  for 


reasons  to  be  discussed  below.  The  current  configuration  is  more  "conservative" 
with  regard  to  the  confinement  of  the  main  flow  gas,  and  the  general  conclusions 
of  the  flow  simulation  experiment  would  therefore  still  be  valid  for  the  present 
design. 

When  the  flow  simulation  experiments  were  performed,  a  number  of  interesting 
results  were  observed.  The  major  result  is  that  with  a  main  flow  Reynolds  number 
greater  than  approximately  100  (for  main : shield : purge  flow  rate  ratios  of  10:5:1) 

the  main  flow  appears  to  be  stable  and  cylindrical  in  form  as  it  passes 
through  the  analysis  region.  Under  these  stable  flow  conditions,  no  recir¬ 
culation  pattern  of  the  main  flow  water  was  observed  in  the  "optics  tubes". 

At  lower  Reynolds  numbers  the  flow  became  unstable  and  main  flow  water  entered 
the  optics  tubes.  In  addition,  when  dye  was  added  to  the  water  of  the  shield 
flow,  a  strong  recirculation  pattern  of  shield  flow  water  was  observed  in 
the  optics  tubes,  and  this  phenomenon  persisted  even  at  the  highest  main  flow 
Reynolds  number  (560)  used  in  the  experiments.  Thus,  it  may  be  concluded  that 
a  simpler  flow  configuration  without  a  shield  flow  would  not  satisfy  the 
analysis  region  design  criteria  since  there  would  be  recirculation  of  absorbing 
gas  in  the  optics  tubes.  The  more  sophisticated  design  is  successful  because 
only  the  non-absorbing  shield  gas  recirculates  in  the  optics  tubes,  while  the 
main  flow  of  absorbing  gas  is  confined  by  the  shield  flow  to  a  cylindrical 
region  in  the  center  of  the  analysis  region.  In  addition,  since  the  shield 
gas  will  nominally  be  at  the  same  temperature  as  the  main  flow,  the  column  of 
absorbing  gas  should  be  essentially  isothermal. 


Although  a  flow  simulation  experiment  of  this  type  can  successfully 
predict  the  bulk  fluid  dynamic  properties  of  an  apparatus,  it  has  a  shortcoming 
in  that  it  can  not  simulate  the  effect  of  diffusion  in  the  gas  phase.  This 
is  important  because  species  in  the  main  flow  which  diffuse  into  the  shield 
flow  could  recirculate  in  the  optics  tubes.  This  was  observed  to  occur  during 
a  series  of  experiments  to  measure  the  temperature  dependence  of  the  spectral 
absorption  coefficients  of  the  band  of  at  which  time  the  apparatus 

was  operated  in  its  original  configuration  with  a  main  flow  nozzle  area 
ratio  of  2.1.  Subsequently,  this  area  ratio  was  increased  to  3.3  (the 
present  design).  In  this  configuration  the  main  flow  was  determined  to 
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be  well  confined  during  measurements  of  the  infrared  band  strength  of  BF. 

In  addition,  it  was  possible  to  verify  visually  the  confinement  of  the  main 
flow  by  observing  in  the  analysis  region  the  afterglow  of  an  air  discharge 
under  the  flow  conditions  used  in  the  BF  experiments.  This  flow  visualization 
experiment  confirms  that  the  apparatus  is  capable  of  successfully  confining 
the  species  in  the  main  flow  as  they  pass  through  the  analysis  region.  However, 
when  different  flow  conditions  are  employed,  it  would  be  advisable  to  repeat 
the  flow  visualization  experiment  for  those  conditions  in  order  to  verify  the 
confinement  of  the  species  in  the  main  flow. 

B.3.  Flow  tube  diagnostics 

In  the  course  of  performing  kinetic  and  spectroscopic  studies,  a  number 
of  key  parameters  which  characterize  the  gases  must  be  accurately  known. 
Calibration  measurements  to  determine  the  temperature,  pressure,  flow  rates, 
and  flow  development  of  the  gases  in  the  flow  reactor  are  described  in  this 


section. 


Radial  Velocity  Distribution  of  the  Flow 

Lair.inar  flow  exists  in  the  main  flow  as  long  as  the  Reynolds  number 

is  below  2000.  This  is  true  for  all  present  applications  of  the  apparatus. 

When  the  carrier  gas  is  introduced  through  the  mullite  arrays,  the  flow  has  a 

uniform  radial  velocity  distribution  (plug  flow).  As  the  gas  flows  down  the 

tube,  its  velocity  profile  eventually  becomes  parabolic  (Poiseuille  flow) . 

In  analyzing  kinetic  data,  the  corrections  for  diffusion  depend  on  the  extent 

to  which  the  flow  has  developed.  Since  the  most  accurate  solutions  to  the  flow 

equations  which  account  for  both  radial  and  axial  diffusion,  as  well  as  wall 
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removal  assume  completely  developed  Poiseuille  flow,  it  is  desirable 

to  operate  in  this  regime.  The  inlet  length  required  to  fully  develop 
parabolic  flow  is1^ 

Z  =  R  d/16  (3) 

e 

where  R  is  the  Reynolds  number  and  d  the  tube  diameter.  In  the  present 
e 

apparatus  the  distance  from  the  laminar  discs  to  the  position  where  reactants 

are  injected  for  kinetic  measurements  is  "-25  cm.  Thus,  fully  developed  flow 

will  present  if  Rg  <  55.  Typically,  Rg  <  10  in  our  kinetics  experiments  so 

that  the  requirements  should  be  satisfied. 

To  verify  this,  we  constructed  a  Pitot  tube,^  which  measures  the  dynarr.i, 

pressure  of  the  gas  as  a  function  of  radial  position.  For  a  compressible  fluid, 

19 

this  pressure  can  be  related  to  the  flow  velocity  by 
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(4; 


where  Pq  is  the  static  pressure,  P  is  the  dynamic  pressure,  Y  the  heat 
capacity  ratio,  and  M  the  gram  molecular  weight.  For  He  at  294  K,  a  Reynrlc.- 
number  of  30,  and  at  a  distance  of  ~25-30  cm  from  the  mullite  arrays,  the  cate 
shown  in  Fig.  4  were  obtained.  The  circles  are  the  measured  velocities  and  the 
solid  line  is  a  parabola  normalized  to  the  center  velocity.  Poiseuille  flow 
has  been  clearly  established. 

Gas  Temperature  Measurements 

In  any  spectral  or  kinetic  measurement,  the  gas  temperature  must  be 

known.  The  problems  encountered  in  measuring  the  gas  temperature  in  a 

variable  temperature  flow  reactor  are  twofold.  First,  the  gas  must  be 

heated  to  the  desired  operating  temperature  before  it  reaches  the  reaction 

and  measurement  region  of  the  tube  and  second,  an  ordinary  thermocouple 

receives  more  energy  via  radiation  from  the  walls  than  from  the  gas,  leading 

to  an  erroneously  high  reading.  The  latter  problem  has  been  overcome  by  using 

3  20 

a  dual  shielded  thermocouple  arrangement  ’  which  was  used  to  calibrate  the 
temperature  controller  readings  as  a  function  of  gas  Reynolds  number.  As  men¬ 
tioned  earlier,  gas  heating  is  accomplished  with  two  mullite  multihole  arrays. 
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Moffat20  has  shown  that  under  laminar  flow  conditions  and  with  a  properlv 
designed  thermocouple,  the  major  source  of  error  in  the  reading  arises  fror 
radiative  heat  transfer  between  the  TC  junction  and  the  reactor  walls.  This 
error  is  reduced  by  heat  shields  at  least  as  much  as  a  factor  of  (n+1),  where 
n  is  the  number  of  concentric  shields  around  the  thermocouple.  By  comparing 
the  readings  from  a  doubly  shielded  and  unshielded  thermocouple,  this  error 
term  can  be  calculated  and  the  true  gas  temperature  measured. 

Such  a  device  was  constructed,  using  W-5%Re/W-26%Re  thermocouples  (,010"dia.), 
referenced  to  0°C,  and  alumina  tubular  shields.  The  inner  shield  is  12  mm  OD  and 
the  outer  one  57  mm  OD.  Both  are  0.5  mm.  thick  and  7.5  cm  long.  Support  is  provide 
by  very  thin  alumina  tubes  used  as  struts  to  minimize  any  conductive  heat  transfer. 
The  unshielded  thermocouple  is  -5  cm  upstream  of  the  shielded  TC. 

The  probe  was  inserted  into  the  flow  tube  at  a  position  midway  through 
the  kinetic  measurements  region.  The  corrected  gas  temperature  was  compared 
with  the  reading  obtained  from  the  thermocouple  on  the  heater  controller 
as  well  as  with  the  "nominal"  value  obtained  by  reading  the  scale  on  its 
temperature  meter.  This  was  done  as  a  function  of  gas  Reynolds  number  fusing 
Ar  and  He  at  various  flow  settings)  for  five  elevated  temperatures.  The 
precision  of  the  TC  readings  was  i  0. 1  K  with  an  accuracy  estimated  at  1  2K. 

The  controller  could  be  read  to  + 5 K.  The  difference  between  the  nominal 
readings  and  the  measure’d  gas  temperatures  are  shown  in  Fig.  5.  For  Re  <  200, 
the  true  gas  temperature  is  ca.  6  K  lower  than  the  nominal  values.  This  varies 
by  +  10  K  for  various  temperatures,  but  shows  no  particular  trend.  Thus  for 
<  200,  the  gas  temperature  is  known  to  t  10  K.  The  difference  between  the 
heater  TC  and  measured  gas  temperatures  exhibits  the  same  behavior. 
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As  Rg  approached  500,  the  difference  between  the  nominal  and  true  values 
precipitously  increased.  This  most  likely  occurs  because  the  smaller  diffu¬ 
sion  and  faster  flow  speeds  at  these  Reynolds  numbers  allowed  insufficient  tim.i. 
for  the  gas  to  be  heated.  The  addition  of  a  third  mullite  array  would  probably 
extend  the  range  of  Rg  where  the  temperature  is  well  characterized.  By  providing 
addition  surface  area  over  which  the  gas  must  pass,  more  heat  can  be  transferred. 

Flow  Rates  and  Gas  Pressures 

Pressures  are  measured  using  an  MKS  Baratron  Model  222AHS100  (1.5T.  accuracy) 
or  MKS  Baratron  Models  310-BHS10,  301-BHS100  (.08/.  accuracy).  The  injector  loo; 
is  used  to  measure  the  static  pressure  inside  the  flow  tube.  This  was  compared 
to  other  probes  in  the  arms  and  entrance  to  the  tube  and  all  readings  agreed 
to  within  1%.  Gas  volumetric  flow  rates  are  obtained  using  Matheson  rotameter;- 
which  are  periodically  calibrated  for  each  gas  used. 


C.  Measurement  Capabilities 
C.l.  Mass  spectrometer 

A  molecular  beam  mass  spectrometer  may  be  interfaced  to  the  high 
temperature  flow  reactor  in  order  to  measure  the  number  density  of  species 
in  the  main  gas  flow.  A  cross  sectional  view  of  the  mass  spectrometer  in 
a  horizontal  plane  is  shown  in  Fig.  6.  A  novel  feature  of  the  design  is 
the  capability  to  measure  the  mass  spectra  of  gases  sampled  from  two  inde¬ 
pendent  apparatus  located  at  right  angles  to  each  other.  Currently,  the 
apparatus  which  is  located  in  the  axial  beam  position  (at  the  left  of  the 
figure)  is  an  atmospheric  pressure,  pulverized  coal  burner,  for  which  the 

mass  spectrometer  measures  the  neutral  and  ionic  species  in  the  burner 
21 

flame.  This  burner  will  not  be  discussed  further  in  the  present  paper. 
However,  it  should  be  noted  that  the  axial  beam  skimmer  and  the  two  inner 
chambers  may  be  electrically  biased  individually  in  order  to  facilitate  the 
measurement  of  ionic  species.  (For  simplicity,  the  additional  ion  optics 
are  not  shown  in  the  figure.) 

The  high  temperature  flow  reactor  is  located  in  the  crossed  beam 

position  (at  the  top  of  the  figure).  From  either  of  the  two  apparatus 

positions,  the  instrument  operates  as  a  conventional  molecular  beam  mass 
22 

spectrometer.  There  are  three  stages  of  differential  pumping.  The  first 
(nozzle-skimmer)  chamber  is  pumped  by  a  16  in.  diffusion  pump  (Varian  type 
HS16),  the  middle  chamber  is  pumped  by  a  4  in.  diffusion  pump  (Varian  type 
VHS-4),  and  the  mass  spectrometer  chamber  is  pumped  by  a  6  in.  diffusion 
pump  (Varian  type  VHS-6) .  A  three  position  gate  valve  mounted  on  the  innermost 
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chamber  permits  operation  of  the  mass  spectrometer  in  one  of  three  modes: 


cross  beam  entrance  and  exit  apertures  open  and  axial  beam  entrance  aperture 
closed  (so  that  the  mass  spectrometer  can  be  kept  evacuated  with  the  outer  two 
chambers  at  atmospheric  pressure).  Choppers  (Bulova  Corp.)  located  on  the  gate 
valve  assembly  permit  the  beams  to  be  modulated  before  they  enter  the  mass  spec¬ 
trometer  chamber,  thereby  permitting  the  use  of  synchronous  detection  techniques 
meter  chamber,  thereby  permitting  the  use  of  synchronous  detection  techniques 
for  discrimination  against  background  gases.  The  mass  spectrometer  proper  is 
a  standard  commercial  unit  (Extranuclear  Laboratories,  Inc.)  except  for  the 
ionizer  (model  041-2)  which  has  been  modified  to  accept  a  cross  beam,  as  well 
as  an  axial  beam. 

The  present  molecular  beam  mass  spectrometer  is  a  recent  addition  to 
the  apparatus,  and,  therefore,  its  ultimate  capabilities  have  yet  to  be 
realized.  At  the  present  time,  when  it  is  operated  with  completely  resolved 
mass  peaks  (Am  *  *s  amu),  it  is  capable  of  detecting  (signal-to-noise  ratio  = 

1,  1  sec  time  constant  on  lock-in  amplifier)  species  at  approximately  4  ppm 
for  the  cross  beam  at  a  source  pressure  of  2  Torr  and  approximately  1  ppm  for 
the  axial  beam  at  a  source  pressure  of  one  atmosphere. 


C.2.  Infrared  absorption 

The  high  temperature  flow  reactor  has  been  designed  with  the  capability 
to  perform  quantitative  measurements  of  infrared  spectral  absorption  coefficients. 
These  measurements  may  be  performed  either  at  low  spectral  resolution  using  a 
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broadband  radiation  source  and  monochromator  or  at  very  high  spectral 
resolution  using  a  tunable  diode  laser  as  the  radiation  source.  The 
infrared  optical  system  is  divided  into  two  major  subassemblies :  one 
Inside  the  flow  reactor  vacuum  jacket  and  one  external  to  the  vacuum 
jacket. 

The  internal  infrared  optical  system  is  shown  in  Fig.  7,  which  is 

a  cross  sectional  view  of  the  flow  reactor  in  the  analysis  region.  The 

internal  optics  are  composed  of  three  spherical  mirrors  (43.0  cm  radius 

23 

of  curvature)  which  form  a  multi-pass  absorption  cell  (White  cell). 

The  White  cell  is  capable  of  producing  approximately  100  passes  through 
the  absorbing  gas  in  the  main  flow,  but  it  is  typically  operated  in  the 
range  of  28-52  passes,  which  is  a  compromise  between  increased  path  length 
(optical  depth)  and  increased  reflectance  losses  as  the  number  of  passes  is 
increased.  The  height  of  the  "envelope"  of  infrared  beams  passing  through 
the  main  flow  is  typically  on  the  order  of  2  cm,  as  compared  with  the  main 
flow  diameter  of  approximately  4  cm. 

The  design  of  the  low  spectral  resolution,  double  beam  external 
optical  system  is  shown  in  Fig.  8.  The  infrared  radiation  source  is  a 
Perkin- Elmer  Corp.  "Opperman  Source"  which  is  used  in  their  current  infrared 
spectrometers.  The  radiation  impinges  on  a  mirror-bladed ,  50  percent  duty 
cycle  chopping  wheel  (Valtec  Corp.)  which  alternately  allows  it  to  pass  into 
the  White  cell  or  reflects  it,  thereby  bypassing  the  White  cell.  In  Fig.  8, 
M2,  M3,  M4,  M7  and  M8  are  flat  mirrors;  Ml  and  M9  are  spherical  mirrors;  and 
M5  and  M6  are  off-axis  paraboloid  mirrors.  The  beam  splitter  is  made  of 
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uncoated  zinc  selenide,  the  monochromator  is  a  Spex  Minimate,  and  the 
detector  is  a  Santa  Barbara  Research  Corp.  HgCdTe  photoconductive  detector. 

In  front  of  this  is  placed  an  order-sorting  (long-wavelength  pass)  filter 
to  avoid  interference  from  higher  order  diffraction  in  the  monochromator. 

Since  the  optical  system  consists  of  reflective  optics  (except  for  the 
zinc  selenide  beam  splitter),  it  is  achromatic,  and  it  can  be  aligned  with 
the  use  of  a  helium-neon  laser. 

The  design  of  the  high  spectral  resolution  external  optical  system  is 
shown  in  Fig.  9.  This  system  shares  many  of  the  components  of  the  low 
spectral  resolution  system.  The  tunable  diode  laser  is  mounted  in  a  closed 
cycle  refrigerator  (both  from  Laser  Analytics,  Inc.).  The  laser  radiation 
is  focused  by  a  KRS5  lens  on  the  monochromator  slit,  and  the  monochromator 
is  used  to  isolate  a  single  laser  mode.  Before  entering  the  White  cell, 
the  laser  radiation  is  chopped  so  that  lock-in  detection  techniques  can 
discriminate  against  blackbody  radiation  from  the  heated  flow  reactor  which 

impinges  on  the  HgCdTe  detector.  With  the  diode  laser,  either  direct  absorp- 
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tion  or  derivative  absorption  techniques  may  be  used  to  measure  the  absorp¬ 
tion  cross  sections  of  individual  spectral  lines.  This  capability  of 
detecting  individual  lines,  as  well  as  the  relatively  high  laser  power, 
permits  an  increase  in  measurement  sensitivity  compared  to  low  spectral 

resolution  detection  techniques.  At  low  spectral  resolution,  a  BF^  number 
13  -3 

density  of  '3  x  10  cm  may  be  detected  using  the  band  (signal-to-noise 
■  1,  time  constant  ■  0.4  sec).^  On  the  other  hand,  at  high  spectral  resolution 
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11  -3 

a  BF  number  density  of  ~10  cm  may  be  detected  by  direct  absorption 

10  -3 

measurements,  and  a  BF  number  density  of  ~2  x  10  cm  may  be  detected 

using  derivative  absorption  techniques  (signal-to-noise  =  1,  time  constant 

_  .  .12 

■  0.4  sec). 

C.3 .  Laser- induced  fluorescence 

The  use  of  laser-induced  fluorescence  extends  the  capability  and 

sensitivity  of  detection  to  species  such  as  NH0.  A  nitrogen-pumped  dye 

laser  (Molectron  DL14P)  is  introduced  through  the  bottom  arm  of  the  flow 

tube  (Fig.  10).  In  this  configuration,  the  resonance  lamp  is  removed  and 

Brewster  windows  added  to  the  arms  along  the  laser  path  to  prevent  beam 

scattering.  The  visible  laser  beam  is  expanded  to  12  ram  diameter  to  avoid 

any  possible  saturation  effects.  (In  the  UV,  this  is  unnecessary  because 

that  laser  power  is  lower.)  After  passing  through  the  flow  tube,  a  portion 

of  the  beam  is  split  off  into  an  ITT  FW-114A  fast  vacuum  photodiode  which 

is  used  as  a  power  monitor.  The  fluorescence  is  imaged  into  a  photomultiplier 

which  is  perpendicular  to  the  path  of  the  laser  beam. 

13  15 

For  each  laser  pulse,  there  are  10  -  10  photons.  Even  with  good 

baffling,  scattered  light  may  create  severe  background  problems.  To  avoid 

this,  a  gated  detection  system  is  used,  and  is  controlled  by  a  master 

25 

electronics  console,  designed  and  built  at  Aerodyne.  This  unit  initiates 

the  laser  pulse,  and  after  the  desired  delay  time,  gates  the  FMT  on  (20  ns 
26  27 

risetime).  '  The  fluorescence  is  amplified,  integrated,  and  fed  to  a 
peak  detector,  where  the  signal  is  held  until  read  into  a  Prime  400  computer 


via  an  IBM  7406  interface.  The  unit  controls  all  delay  and  data  integration 

times  in  increments  of  10  ns.  A  parallel  channel  repeats  this  process  for 

the  laser  power  monitor.  The  laser  is  fired  on  alternate  clock  periods  and 

the  signal  (difference  between  laser-on  and  laser-off  values,  normalized  to 

the  laser  power  for  each  pulse)  is  computed  and  displayed  in  real  time  on  a 

9 

video  terminal.  For  a  typical  species  such  as  OH,  a  concentration  of  ~10 

_3 

molecules  cm  can  be  detected  in  10  seconds  of  counting  with  signal/background 
=  1. 

C.4.  Resonance  fluorescence 

A  variety  of  atomic  and  molecular  species  can  be  detected  using  resonance 
fluorescence  techniques.  A  general  schematic  of  this  configuration  is  illus¬ 
trated  in  Figs.  2  and  10.  The  quartz  resonance  lamp  is  made  by  flowing  the 
desired  species,  entrained  in  an  inert  carrier,  through  a  2450  MHz  microwave 
discharge  and  focussing  the  resulting  emission  into  the  flow  tube. 

Fluorescence  is  imaged  from  the  central  12  mm  diameter  region  of  the 
analysis  region  of  the  flow  tube  onto  a  photomultiplier  configured  for  photon 
counting.  The  PMT  is  mounted  on  an  arm  perpendicular  to  the  resonance  lamp. 

The  counts  are  amplified  and  discriminated  by  a  Pacific  Precision  Instruments 
Model  3262/AD6  amplifier-discriminator  and  fed  to  a  Monsato  1505A  Counter-Timer. 
Appropriate  interference  filters  are  used  to  isolate  the  emission  lines  from 
unwanted  scattered  light,  which  can  be  very  intense  at  elevated  temperatures. 

When  detecting  oxygen  or  hydrogen  atoms,  a  solar  blind  EMI  G-26E315 
photomultiplier  is  used,  with  optics  constructed  from  MgF2*  For  all  other 
species  (BF,  NO,  OH,  NO2  etc.),  quartz  optics  and  a  Hammamatsu  R763P  phototube 
(160-830  nm  response)  are  used. 
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As  shown  in  Figs.  2  and  10,  the  optical  paths  are  protected  from  stray 
light  by  black  anodized  aluminum  baffles.  Opposite  the  PMT  is  a  blackened 


light  trap,  which  suppresses  scattering  light  so  that  the  PMT  views  no 

direct  source  of  emission  other  than  the  fluorescing  gases.  It  is  estimated 

1 2 

that  scattered  light  is  attenuated  by  ~10  .  For  a  typical  species  such 

10  -3 

as  NO,  a  concentration  of  — 1.3  x  10  molecules  cm  can  be  detected  in 
10  seconds  of  counting  with  a  signal/background  =  1. 

C.5.  Ultraviolet  absorption 

Species  which  have  moderate  absorption  coefficients  in  the  UV  may  be 
detected  by  direct  absorption  of  continuum  radiation  from  a  xenon  arc  lamp 
(Varian  Model  VIX-500  UV) .  The  lamp  has  a  total  output  of  12  W  integrated 
from  200  to  400  nm.  Although  the  intensity  decreases  rapidly  below  220  nm, 
there  is  still  a  useful  continuum  from  190  to  200  nm  of  1  mW/nm.  As  shown 
in  Fig.  7,  the  light  is  collimated  by  an  internal  parabolic  reflector,  passes 
through  the  center  of  the  detection  region,  and  is  focused  by  a  quartz  lens 

onto  the  slit  of  a  1  m  monochromator  (McPherson  Model  2051).  The  1200  line/mir. 

o  o 

grating  is  blazed  at  3000  A  and  gives  a  resolution  in  first  order  of  0.1  A 

with  10  pm  slits.  A  beam  splitter  deflects  a  portion  of  the  lamp  output  onto 
a  reference  photomultiplier  tube  so  that  intensity  fluctuations  may  be  com¬ 
pensated  by  ratioing  the  two  beams.  The  optical  axis  for  the  UV  absorption  is 
perpendicular  to  the  infrared  axis  permitting  simultaneous  measurements  by 

the  two  techniques.  For  a  typical  species  such  as  BF,  a  number  density  of 
11  -3 

—2  x  10  cm  can  be  detected  with  a  signal-to-noise  *  1  with  a  time  constant 
of  10  sec.  This  sensitivity  is  largely  limited  by  fluctuations  in  the  arc  lamp. 
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C.6.  Chemiluminescence 


Visible  and  ultraviolet  chemiluminescence  spectra  are  taken  by  imaging 
emission  from  the  central  region  of  the  flow  axis  onto  the  entrance  slit  of 
the  monochromator  using  the  same  optical  arrangement  as  in  Fig.  7.  Spectra 
are  recorded  with  the  IBM  7406  Interface  and  Prime  400  Computer  using  either 
analog  or  digital  data  from  the  phototube  while  wavelength  is  scanned  by  a 
stepping  motor  drive. 
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II .  OPERATION 


The  high  temperature  flow  reactor  has  been  used  for  several  spectroscopic 
and  kinetic  studies.  These  studies  will  be  summarized  briefly  in  order  to 
illustrate  the  versatility  of  the  apparatus. 

The  temperature  dependence  of  the  infrared  spectral  absorption  coefficients 
of  the  l> band  of  BF^  was  measured  over  the  temperature  range  of  298  to  1440  K.* 
The  BF^  was  simply  mixed  with  the  carrier  gas  in  the  main  flow,  and  the  infrared 
absorption  coefficients  were  measured  at  low  spectral  resolution.  The  BF^ 
number  density  in  the  flow  tube  was  obtained  from  the  gas  flow  rates  and  the 
pressure  in  the  analysis  region. 

-2 

The  infrared  band  strength  of  BF  was  determined  to  be  630  l  200  cm 

amagat  1  by  measuring  the  integrated  line  strength  of  a  single  vibration- 

12 

rotation  line  at  high  spectral  resolution.  The  BF  was  produced  by  passing 
a  mixture  of  BF^  and  He  carrier  gas  through  a  microwave  discharge.  The  BF 
number  density  in  the  analysis  region  was  obtained  by  measuring  its  absorption 
of  ultraviolet  radiation  along  a  path  perpendicular  to  the  path  of  the  reaction 
of  the  infrared  radiation. 

The  rate  constant  for  the  reaction  of  OH  +  NH^  was  measured  over  the 

26  —12 
temperature  range  of  294  to  1075  K  and  found  to  be  (5.41  ±  0.86)  x  10 

-1  3  -1  -1 

exp (-2120  ±  143  cal  mole  /RT)  cm  molecule  sec  .  The  OH  was  produced  by 
reaction  of  H  (generated  in  a  microwave  discharge)  with  and  the./ NH^  was 

introduced  into  flow  downstream  through  the  ring  injector.  The  progress  of 
the  reaction  was  monitored  in  the  analysis  region  by  measuring  the  OH  number 
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density  by  resonance  fluorescence  or  by  measuring  the  N^  product  number 
density  by  laser  induced  fluorescence. 

Finally,  the  rate  constant  for  the  reaction  of  NH^  +  NO  was  measured 

29 

over  the  temperature  range  of  294  to  1200  K.  A  least  squares  fit  of  the 

-4  -2  46 

experimental  data  yielded  the  expression  k(T)  =  (1.95  x  10  )T  '  exp(-938/T) 

3  -1-1 

cm  molecule  sec  .  The  NH^  was  produced  by  reaction  of  F  (generated  in  a 
microwave  discharge)  with  NH^,  and  the  NO  was  introduced  through  the  ring 
injector.  In  the  analysis  region  the  NH2  reactant  number  density  was  measured 
by  laser  induced  fluroescence,  the  OH  product  number  density  was  measured  by 
laser  induced  fluorescence  and  resonance  fluorescence,  and  the  number  densities 
of  all  the  species  in  the  main  flow  were  measured  using  the  molecular  beam 
mass  spectrometer. 
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FIGURE  CAPTIONS 


Figure 

1.  Block  diagram  of  major  apparatus  components. 

2.  Cross  sectional  view  of  high  temperature  flow  reactor. 

A  -  alumina  flow  tube,  B  -  water  cooled  vacuum  jacket, 

C  -  refractory  insulation,  D  -  ring  injector,  E  -  six  arm 
cross  piece  to  support  flow  tubes,  F  -  discharge  lamp, 

G  -  microwave  discharge  cavity,  H  -  quartz  lens,  I  -  flow 
of  gas  to  pumps,  J  -  optical  tube  perpendicular  to  drawing, 

K  -  mass  spectrometer  sampling  orifice,  L  -  optics  tube 

purge  gas  (1  of  4),  M  -  kanthal  heating  elements,  N  -  multihole 

mullite  array  (heat  exchanger),  0  -  shield  gas,  P  -  carrier  gas, 

Q  -  reactant  gases. 

3.  Cross  sectional  view  of  analysis  region. 

4.  Measured  velocity  profile  in  flow  tube.  Solid  line  is  parabola 
normalized  to  data  point  on  centerline. 

5.  Difference  between  nominal  readings  and  measured  gas  temperatures. 

6.  Cross  sectional  view  of  molecular  beam  mass  spectrometer. 

A  -  nozzle-skimmer  chamber,  B  -  middle  chamber,  C  -  mass  spectro¬ 
meter  chamber,  D  -  sampling  nozzles,  E  -  skimmers,  F  -  choppers. 

G  -  gate  valve  assembly,  H  -  mass  spectrometer,  I  -  nylon  flanges 
permit  electrical  biasing  for  detection  of  ions  from  coal  burner. 

7.  Cross  sectional  view  of  analysis  region  with  infared  absorption 
and  ultraviolet  absorption  diagnostics.  Elements  external  to 
vacuum  jacket  not  shown  to  scale.  A  -  vacuum  jacket,  B  -  six 
arm  cross  piece,  C  -  white  cell  mirrors,  D  -  windows,  E  -  arc 

lamp,  F  -  beam  splitter,  G  -  interference  filters,  H  -  photomultipliers 
I  -  mirror,  J  -  quartz  lens,  K  -  chopper,  L  -  monochromator. 

8.  View  of  low  spectral  resolution  infrared  optics.  A  -  source, 

B  -  chopper,  C  -  beam  splitter,  D  -  monochromator,  E  -  detector, 

M1-M9  -  mirrors. 

9.  View  of  high  spectral  resolution  infared  optics.  A  -  tunable  diode 
laser,  B  -  lens,  C  -  monochromator,  D  -  chopper,  E  -  detector, 

M1-M8  -  mirrors. 

10.  Cross  section  view  of  analysis  region  with  laser  induced  fluorescence 
diagnostics.  A  -  vacuum  jacket,  B  -  six  arm  cross  piece,  C  -  pulsed 
dye  laser,  D  -  beam  expander,  E  -  mirrors,  F  -  Brewster  angle  windows, 

G  -  baffle  assemblies,  H  -  photodiode,  I  -  light  trap,  J  -  quartz 
lens,  K  -  interference  filter,  L  -  photomultiplier. 
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3.  MEASURED  AND  PREDICTED  BAND  MODEL  PARAMETERS 
OF  BF3(v3)  AT  HIGH  TEMPERATURES 

by 

Michael  E.  Gersh  and  Lawrence  S.  Bernstein, 
and  Lauren  M.  Peterson 

ABSTRACT 

Two  experimental  techniques  were  used  to  determine  the  high  temperature 
infrared  absorption  coefficients  of  the  v3  band  of  BF^ .  A  high  temperature 
flow  tube  fitted  with  a  multipass  infrared  absorption  cell  was  used  to 
determine  the  BF^  absorption  coefficients  at  298,  585,  870,  1150,  and  1440  K. 
BF^  was  also  used  in  a  flame  emission-absorption  apparatus  to  determine  the 
BF3  absorption  coefficients  at  2400  K. 

A  generalized  formulation  for  predicting  the  S/d  and  1/d  band  model 
parameters  of  symmetric  top  perpendicular  bands  of  MX3  molecules  is  presented. 
The  effect  of  including  nuclear  spin  statistics  in  computing  the  S/d  and  1 /d 
band  model  parameters  is  considered.  It  is  shown  that  nuclear  spin  statistics 
effect  only  the  1/d  parameter  and  are  significant  only  for  the  case  of  zero 
nuclear  spin.  A  number  of  practical  considerations  which  arise  in  designing 
a  computer  algorithm  to  calculate  the  S/d  and  1/d  parameters,  are  discussed. 

Comparisons  of  the  theoretical  and  experimental  BF3  absorption  coefficients 
are  made.  Good  agreement  between  theory  and  experiment  at  all  temperatures 
was  obtained  by  adjusting  only  two  vibrational  energy  coupling  constants, 

X23  ®nd  Predicted  S/d  and  1/d  band  model  parameters  at  300,  600,  900, 

1200,  1500,  2000,  2500,  and  3000  K  are  presented. 
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I .  INTRODUCTION 


For  diatomic  and  linear  triatomic  molecules,  theoretical  predictions  and 

experimental  determinations  of  high  temperature  infrared  absorption  coeffi- 
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cients  are  in  close  agreement.  Our  present  knowledge  of  high  temperature 

absorption  coefficients  for  polyatomic  molecules  is  less  complete.  As  a  first 
venture  into  high  temperature  polyatomic  molecular  spectroscopy  we  have  per¬ 
formed  a  combined  theoretical  and  experimental  study  of  the  band  of  BF^. 

The  vibrational  and  rotational  constants  of  the  ground  and  first  excited 

vibrational  states  of  BF^  have  been  experimentally  determined  by  previous 
14-35 

efforts.  In  addition,  the  general  theoretical  description  of  fundamental 

vibrational-rotational  band  contours  for  symmetric  top  molecules  is  well 

established"*6-^1  and  has  been  shown  to  adequately  describe  the  room  temperature 

2A37  14  15  16  1G  0 

infrared  spectrum  of  BF^  at  both  low  ’  and  high  spectral  resolution.  ’  ’  ’ 

This  work  provides  the  foundation  for  the  present  study. 

Two  complementary  techniques  were  used  for  obtaining  high  temperature 
infrared  absorption  coefficients.  These  are:  (1)  absorption  measurements  in 
a  high  temperature  flow  reactor  and  (2)  flame  emission-absorption  measurements. 

In  the  flow  reactor,  gaseous  mixtures  are  conductively  heated  to  temperatures 
up  to  1500  K,  and  an  aerodynamically  confined  flow  is  established  which  is 
Intersected  by  a  vindowless  multipass  infrared  absorption  cell.  In  the  flame 
studies,  BF^  emission-absorption  measurements  were  performed  downstream  of 
the  burner  exit.  A  temperature  of  2400  K  was  attained  in  the  flame  measure¬ 
ments.  The  experimental  data  provide  a  complete  picture  of  the  temperature 
dependent  band  shape  of  the  band  of  BF^  from  room  temperature  to  2400  K. 
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These  data  were  used  to  validate  a  generalized  band  model  for  symmetric  top 
perpendicular  bands  of  MX^  molecules,  which  is  applicable  to  both  the  planar 
and  pyramidal  geometries. 
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II.  GLOSSARY  OF  NOTATION 


JK 

<A  > 
v 


<  B  > 

v 


c 

C, 


absorptivity 

equilibrium  rotational  constant 
Honl-London  intensity  factor 
average  rotational  constant 

rotational  constant  for  v'th  vibrational  state 

equilibrium  rotational  constant 

rotational  constant  for  v'th  vibrational  state 

average  rotational  constant 

speed  of  light 

second  radiation  constant 


D 


J 


D 


JK 


D 


K 


P,  R 

d  (1) 


centrifugal  distortion  constants 


line  spacing  for  P,  R-branch  transitions 


line  spacing  for  Q-branch  transitions 


E(v,  1,  J,  K)  total  vibration-rotation  energy 


E  Baximuir.  allowed  vibrational  energv  for  computational  convergence 

Bax 


positive  vibrational  energy  term  for  computational  convergenc 


vibrational  angular  momentum  energy  constants 

nuclear  statistical  weight  for  vr'th  transition 
Planck's  constant 
nuclear  spin 

unattenuated  and  attenuated  photon  beams  respectively 

rotational  quantum  number 
rotational  quantum  number 
Boltzman's  constant 

spectral  absorption  coefficient 

physical  absorption  path  length 

vibrational  angular  momentum  quantum  numbers 

spectral  radiance 
rotational  quantum  number 

rotational  quantum  number 

pressure 
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sj“> 

P,  R 
S(w) 
v 

S/d (f) 


vr 


<S/d> 

I 


u 


total  Integrated  band  strength 

total  Integrated  band  strength  for  mode 

Integrated  Q-branch  line  strength 

integrated  P,  R-branch  line  strength 
spectral  absorption  coefficient 

line  strength  of  r'th  transition 

band  averaged  absorption  coefficient 
temperature 

optical  path 


W  (w,  T)  Planck  blackbody  function 
Uc 


z 


v 


Z 

r 


} 


l/d(w) 

<l/d> 


<l/d>' 


vibrational  and  rotational  partition  functions 

line  density  function 

band  averaged  line  density  including  nuclear  spin  statistics 
band  averaged  line  density  not  Including  nuclear  spin  statistics 
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Creek  Symbols 


vibration-rotation  energy  coupling  terms 


y 


!3 


C4 


} 


V 


vi 


•\ 


V 


■> 


\ 


3 


V 


4 


J 


spectral  line  vidth 


Coriolis  coupling  terms 


v'th  vibrational  state 


vibrational  quantum  numbers  for  four  modes  of  an  M>.^  sywi.ie t r 


c  rotational  symmetry  number 

t  transmittance 


if  tCT 


xir  x22’  1 


x33*  X44* 


Xj2»  Xj3>  v  constants  for  quadratic  vibrational  energy  terms 


xu.  x23. 


x24’  x34  J 
X1 1 1 »  x222’ 


x333*  X444 


constants  for  cubic  vibrational  energy  terms 


‘3’  “4  J 


wavenumber 


harmonic  vibrational  energy  constants 


P,  R-branch  transition  energy 
Q-branch  transition  energy 

purely  vibrational  contribution  to  origin  of  v'th  band 

origin  of  Q-branch  for  v’th  band 

modification  of  u.  , 
vJ 

origin  of  P,  R-branches  for  v'th  band 

modification  of  u  ,, 
vh 


Au 


difference  of  vibrational  band  center  location  between  v '  and 
v’ th  bands 

effective  band  width 
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III.  EXPERIMENTAL 


The  measurements  were  performed  in  two  different  types  of  apparatus. 

The  first  is  a  high  temperature  flow  reactor  into  which  the  BF^  was  introduced. 
This  was  used  to  measure  the  absorption  spectra  at  temperatures  between  298  K 
and  1440  K.  The  second  is  a  radiative  transfer  apparatus  that  was  used  to 
measure  the  emission  and  absorption  spectra  of  BF^  in  a  2400  K  flame . 

A.  High  Temperature  Flow  Reactor 

A  cross  sectional  view  of  the  Aerodyne  Research,  Inc.  heated  flow 
reactor  is  shown  in  Fig.  1.  In  general,  the  infrared  active  species  of  interest 
are  produced  by  reaction  in,  or  introduced  into,  a  variable  temperature  flow 
reactor  which  can  attain  temperatures  of  1500  K.  Then,  the  flow  is  passe, 
through  an  analysis  region  in  which  a  multi-pass  absorption  cell  (White  Ct 1 1  '> 
is  mounted  perpendicular  to  the  flow  direction.  Final lv,  a  portion  of  the  flc. 
is  sampled  by  a  molecular  beam  mass  spectrometer  before  beinr  exhausted  tc  z. 
large  mechanical  vacuum  pump.  In  the  present  study  the  BF^  number  densitv 
was  obtained  from  total  pressure,  temperature  and  gas  flow  rates.  For  reactive- 
species  produced  in  the  flow  reactor,  their  number  density  can  be  measured 
in  situ  by  the  mass  spectrometer  or  by  laser-induced  or  resonance  fluorescence 
(using  mutually  perpendicular  optics  tubes).  Thus,  one  can  obtain  absolute 
values  for  the  absorption  coefficients  without  assuming  that  the  reactive 
species  in  the  flow  tube  have  achieved  local  thermodynamic  equilibriur  . 

The  detailed  design  of  the  flow  reactor  is  discussed  in  Ref.  43. 
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The  design  of  the  gas  flow  system  within  the  variable  temperature  flow 
reactor  is  fairly  complicated  due  to  the  constraints  imposed  by  the  necessity 
of  making  infrared  absorption  measurements.  The  crucial  element  of  the  apparatus 
design  is  the  configuration  of  the  spectroscopic  analysis  region.  In  order  for 
one  to  be  able  to  perform  quantitative  absorption  (or  emission)  measurements 
on  the  flow  tube  gases,  three  criteria  must  be  satisfied;  (1)  the  geometry  of 
the  absorbing  gas  flow  must  be  well  defined  in  the  region  where  it  traverses  the 
White  cell;  (2)  the  flow  in  this  region  must  be  stable;  and  (3)  the  temperature  dis¬ 
tribution  of  the  absorbing  gas  flow  must  be  well  defined  (preferably,  the  gas  shoul 
be  isothermal).  To  satisfy  these  criteria,  the  flow  of  infrared  active  gas 
should  ideally  be  confined  in  a  cylindrical  geometry  and  be  located  in  a  stable 
position,  without  any  recirculation  of  absorbing  gas  into  the  optics  tubes 
(i.e.,  the  flow  should  behave  as  if  it  were  confined  by  a  cylindrical  tube). 

For  an  ambient  temperature  flow  with  single-pass  optics,  one  could  simply 
place  optical  windows  on  the  flow  tube;  however,  the  present  and  future  experi¬ 
ments  require  high  absorbing  gas  temperatures  and  a  multi-pass  absorption  cell 
(for  additional  sensitivity).  This  necessitates  the  use  of  a  wall-less  con¬ 
finement  of  the  absorbing  gases  as  they  pass  through  the  analysis  region.  This 
can  only  be  accomplished  by  fluid  dynamic  confinement  of  the  absorbing  gas 
flow. 

The  analysis  region  design  shown  in  Fig.  2  was  considered  to  be  most 
promising.  The  main  flow  containing  argon  carrier  gas  and  BF^  is  heated  as 
It  passes  through  two  2.5  cm  thick  mullite  "honeycomb"  hole  arrays,  and  the 
outer,  coaxial  shield  gas  flow  of  argon  is  heated  as  it  passes  over  the  heater 


elements.  Then,  the  main  and  shield  gas  flows  are  both  passed  through 
concentric  subsonic  nozzles,  after  which  they  join  and  pass  through  the 
analysis  region.  The  area  ratio  of  the  main  flow  nozzle  to  the  shield 
flow  nozzle  is  2:1,  and  the  corresponding  flow  rates  are  also  designed  to 
be  in  a  2:1  ratio  so  that  the  velocities  of  the  two  flow’s  will  be  the  same 
when  they  meet  downstream  of  the  nozzles.  In  addition,  it  was  estimated 
that  the  flow  of  purge  gas  into  the  optics  tubes  should  equal  approximately 
10  percent  of  the  main  flow  and  should  be  evenly  divided  between  the  four 
optics  tubes. 

This  flow  geometry  is  considerably  more  complicated  than  that  used 

44  45 

previously  for  flow  reactors,  ’  and  the  design  criteria  for  the  conf  inener.t 

of  the  absorbing  gas  flow  are  quite  strict.  Thus,  it  was  concluded  that  this 

crucial  element  of  the  flow  system  design  should  be  tested  in  a  flow  simulation 
43 

experiment.  The  result  of  this  experiment  was  that  the  main  flow  should  be 
well  confined  for  a  Reynolds  number  of  greater  than  approximately  100  in  the 
main  flow  tube  and  with  the  flow  rate  ratios  given  above  (main: shield :purge =  1 0 : 5 : 1 ) . 

The  design  of  the  present  external  optical  system  for  the  infrared  absorp¬ 
tion  measurements  is  shown  in  Fig.  3.  This  double  beam  optical  system  has  beer, 
designed  for  low  spectral  resolution  measurements.  The  infrared  radiation 
source  is  a  Perkin-Elmer  Corporation  "Opperman  Source"  which  is  a  mixed  oxide  tube 
heated  by  an  internal  coil  of  resistance  wire.  The  radiation  impinges  on  a 
mirror-bladed,  50  percent  duty  cycle  chopping  wheel  which  alternately  allows 
it  to  pass  into  the  White  cell  or  reflects  it,  thereby  bypassing  the  White  cell. 

M2,  M3,  M4,  M7  and  M9  are  flat  mirrors;  Ml  and  M8  are  spherical  mirrors,  and 
M5  and  M6  are  off-axis  paraboloid  mirrors.  The  beam  splitter  is  made  of  uncoated 
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zinc  selenide,  the  monocbromater  is  a  Spex  Minimate,  and  the  detector  is 
a  Santa  Barbara  Research  Corporation  HgCdTe  photoconductive  detector,  in  front 
of  which  is  placed  an  order-sorting  (long-wavelength  pass)  filter. 

In  general,  the  data  are  acquired  and  analyzed  in  the  following  manner. 

At  a  given  wavelength  the  double  beam  optical  system  is  used  to  measure  the 
difference  in  signal  (Io  -  I)  due  to  the  presence  of  BF^  in  the  gas  flow. 

Then,  the  intensity  of  the  single  beam  passing  through  the  gas  flow  I  is 
measured  in  the  absence  of  BF^.  The  ratio  of  these  measurements  is  related 
to  the  spectral  absorption  coefficient  k  by  Eq.  (5). 

(1  -  I) /I  *  1  -  exp(-ku)  C'l 

o  o 

where  k  has  the  units  cm  1  amagat  1  and  u  is  the  optical  path  which  is  related 
to  the  physical  absorption  path  length  (cm)  by 

u  ■=  £  (P/760)  (273/T)  (,2s* 

In  this  equation  P(torr)  is  the  partial  pressure  of  BF^,  and  T(K)  is  the  temperature 
of  the  flowing  gas.  Thus,  if  i,  P  and  T  are  known,  one  can  use  the  data  to  cal¬ 
culate  the  spectral  absorption  coefficient  k. 

The  measurements  were  performed  under  the  following  conditions.  The 
spectral  scans  covered  the  region  from  6  to  8  pm,  with  a  spectral  resolution 
of  0.02  um.  The  wavelength  scale  was  calibrated  by  comparing  a  water  vapor 


absorption  spectrum  with  the  absorption  peaks  measured  by  Oetjen,  et.  al. 


The  temperatures,  total  pressures,  densities  of  BF^,  and  main  flow  Reynolds 
numbers  are  given  in  Table  I.  The  number  of  White  Cell  passes  was  52,  with 
a  nominal  pathlength  i  ■  4.8  cm.  Finally,  the  BF^  was  obtained  from  Matheson 
and  had  a  stated  purity  of  99.5  percent. 
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B.  Flame  Emission-Absorption  Apparatus 


The  radiative  transfer  apparatus  used  at  the  Environmental  Institute  of 
Michigan  to  obtain  IR  flame  spectra  of  BF^  is  shown  schematically  in  Fig.  4. 

It  is  a  dual  beam,  ratio  recording  instrument  utilizing  a  0.27  m  Perkin-Elner 
Model  99G  spectrometer.  Both  emission  and  absorption  measurement  can  be  made, 
and,  since  the  emission  measurements  can  be  calibrated  using  the  known  radiant 
of  the  1400  K  blackbody  source,  absolute  radiance  measurements  can  be  made  ir 
order  to  determine  the  temperature  of  the  flame.  The  minimum  resolution  of  th 
instrument  was  about  2  cm  \  but  the  slits  were  usually  opened  to  provide  a 
resolution  of  about  5  cm  Using  a  liquid  ni trogen-cooled  PbSnTe  d elector  an 
an  integration  time  of  0.3  sec,  a  signal-to-ncise  ratio  of  about  200  was  obtai 
A  seven  segment  23  cm  water-cooled  multiple  tube  diffusion  burner  was  use 
Each  segment  was  3.3  cm  long  and  had  independent  fuel  ports.  A  segment  consis 
of  a  bundle  of  tubes  through  which  the  fuel  was  introduced  to  the  flare  and 
around  which  the  oxidizer  passed,  such  that  the  mixing  region  was  formed  a:  th 
upper  surface  of  the  tubes.  By  introducing  BF^  into  any  of  the  central  five- 
segments,  isothermal  path  lengths  from  3.3  to  16.3  cm  were  obtained  for  the 
active  species  (BF^) .  It  is  worth  noting  that  a  non-isothermal  condition 
exists  between  the  shielding  flame  (end  segments)  and  the  ambient  air,  but, 
since  there  is  no  BF^  present,  this  is 
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of  no  consequence  to  the  observed  spectrum.  Also,  since  the  mole  fraction  o: 

BF^  was  typically  only  0.06  in  the  active  segments  and  since  BF^  contributes 
only  minimally  in  the  combustion  reaction  to  generate  BOF,  isothermal  conditions 
were  present  at  the  active/shield  flame  interface. 

The  burner  position  was  adjusted  so  that  the  beam  of  blackbodv  radiation 
was  centered  in  the  flame  and  was  located  at  a  height  above  the  burner  fact 
which  was  above  the  visible  reaction  zone,  but  low  enough  so  that  the  mixing 
layer  with  the  ambient  air  had  not  penetrated  a  substantial  distance  into  the 
burner  plume.  Heights  of  about  2  cm  were  used. 

During  a  given  spectral  run,  five  complete  spectral  scans  were  performed 
in  the  following  order: 

1.  An  emission  background  scan  with  no  active 
species  (BF^)  present  in  the  flame 

2.  An  absorption  background  scan  with  no 
active  species. 

3.  An  absorption  scan  with  the  active  species 
present  in  the  central  segments  (1,  3  or  5) 
of  the  7-segment  burner . 

4.  An  emission  scan  with  the  active  species  present. 

5.  An  absorption  scan  using  the  emission  chopper. 


3-16 


The  net  absorption  A  was  obtained  by  ratioing  the  transmissions  t : 


TBX 


(3) 


where  X  denotes  the  species  of  interest,  B  denotes  background,  and  BX 
denotes  the  combined  spectrum  due  to  active  species  plus  background.  In 
order  to  obtain  absolute  spectral  radiance,  L (co >  values,  the  recorded  data 
(voltage  values)  are  compared  to  the  reference  radiance  provided  by  the  1400  K 
blackbody.  Flame  temperatures  were  determined  by  dividing  by  ar.;: 

calculating  the  corresponding  temperature  for  an  equivalent  blackbodv,  W£.  (_  ,  T'l. 


LUO/Aj.Cu)  =  VBE(.,  T) 

As  an  independent  check  on  the  temperature  determined  by  emission/absorr:  ior: 
pyrometry,  the  sodium  D  line  reversal  technique^*®  was  used.  Both  methods 
gave  essentially  identical  temperatures  for  the  BF^  measurements. 

In  order  to  quantify  the  IR  emission  and  absorption  measurements,  it  is 
necessary  to  know  the  density  of  active  species  present  i:.  the  flame.  Assuming 
the  flame  is  in  local  thermodynamic  equilibrium,  chemical  equilibrium  calcula¬ 
tions  can  provide  an  accurate  determination  of  the  quantity  of  reaction  products. 
Given  the  flow  conditions  of  the  reactants  and  values  for  their  heats  of 
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formation,  the  flame  temperature  and  the  mole  fraction  of  reaction  products 

49 

present  in  the  flame  may  be  calculated.  Unfortunately,  this  calculated  value 
of  the  flame  temperature  is  usually  less  reliable  than  experimentally  determined 
temperatures,  since  the  calculation  is  ideal  and  does  not  account  for  any  loss 
mechanisms  which  may  be  present  in  the  combustion  process.  For  this  reason, 
the  flame  temperature  determined  by  emission/absorption  pyrometrv  was  used  as 
an  input  condition  to  a  chemical  equilibrium  program,  and  the  computed  molt 
fractions  of  reaction  products  were  used  in  calculating  the  absorption  coefficient 
and  band  strengths  for  the  BF^.  Once  the  temperature  and  composition  of  the  flart 
are  known  quantitatively,  the  optical  depth,  u,  may  be  determined  in  the  sar t 
manner  as  discussed  for  the  flow  reactor  (see  Eq .  (2)).  The  spectral  absorption 
coefficient  is  determined  from  u  and  the  measurement  transmission,  t  (Ec.  (3)',  r; 


k 


In 

u 


(O' 


The  mole  fraction  of  BF^  in  all  the  data  runs  was  approximately  0.06. 
BOF  and  BO^  were  also  present  in  the  flame,  but  their  concentrations  were 
approximately  one  and  four  orders  of  magnitude  lower  than  that  of  BF^, 
respectively.  The  computed  equilibrium  mole  fractions  of  BOF  and  BO^  as  a 
function  of  input  BF^  mole  fraction  can  be  found  in  Fig.  5.  Since  BOF  is 
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not  expected  to  have  any  bands  in  this  spectral  region  and  since  BC>2  is  so 
much  more  dilute  than  the  BF^,  we  may  ascribe  our  measurements  to  be  solely 
due  to  BF^. 

The  emission/absorption  pyrometry  temperature  measurements  yielded  the 
same  temperature  throughout  the  entire  spectral  band.  This  indicates  a  temp¬ 
orally  stable  and  spatially  isothermal  active  flame.  Examination  of  the  active 
flame  with  BF^  present  revealed  a  light  green  plume  over  the  majority  of  the 
active  burner  segments  and  a  bright  green  intermediate  2one  between  the  active 
and  shielding  flame.  The  bright  green  flame  is  caused  by  B02  emission  and  is 
an  indication  of  a  decreased  amount  of  the  reactant  BF^  due  to  the  combustion 
process.  It  is  reasonable  to  assume  that  some  of  the  BF^  is  transported  fror 
the  active  region  of  the  flame  into  the  shield  due  to  mixing  at  the  active 'shield 
flame  interface,  but  the  net  absorption  and  emission  of  the  flame  should  be  the 
same  if  the  active  flame  is  isothermal  and  the  total  number  density  of  BF^ 
remains  unchanged.  As  is  indicated  by  the  optical  temperature  determination, 
the  flame  is  indeed  isothermal,  but  examination  of  the  chemical  equilibrium 
computations  as  a  function  of  input  BF^  concentration  shows  that  the  total 
number  density  of  BF^  in  the  flame  can  change.  As  the  input  mole  fraction  of 
BF^  is  decreased  (i.e.,  by  transport  into  the  shield),  the  mole  fraction  of 
BF^  in  the  combustion  products  does  not  simply  decrease  in  the  same  amount, 
but  the  reaction  mixture  in  the  shield/active  flame  mixing  region  is  modified 
to  drive  more  of  the  BF^  into  the  other  boron  compounds.  This  process  should 
appear  as  a  nonlinear  deviation  from  Beer's  law  behavior  because  the  decrease 
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in  BF^  optical  depth  due  to  mixing  is  proportional  to  the  concentration  of 
BF^  in  the  active  region.  This  nonlinearity  was  observed  and  corrected  for 
by  determining  a  smaller  effective  length  of  the  active  region.  Subtracting 
1  cm  from  all  active  burner  lengths  gave  the  desired  linearity. 
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IV.  THEORETICAL  PREDICTIONS  OF  BAND  MODEL  PARAMETERS 

In  this  section  general  equations  for  the  S/d  and  1/d  band  model  parameters 
are  derived.  The  derivation  (sections  IV  A-C)  is  for  perpendicular  bands  of 
MX^  symmetric  tops  and  applies  to  both  the  planar  (D^)  and  pyramidal  (C^v) 
geometries.  The  effect  of  nuclear  spin  statistics  on  the  calculation  of  the 
band  model  parameters  is  discussed  in  section  IV  D.  In  section  IV  E  approxi¬ 
mations  to  the  general  expression  for  the  band  model  parameters  are  introduced 
which  significantly  decrease  the  computer  time  required  to  evaluate  these 
quantities  at  high  temperatures. 

A.  Average  Line  Spacing 

The  vibration-rotation  energy  of  an  MX^  symmetric  top  molecule  is 
approximated  by^ 

E(v,  1,  J.  K)  -  tojCVj  4  1/2)  +  u.2(v2  4  1/2)  ♦  w3(v3  *  1) 

4  4  i>  4  Xji(vj  4  1/2>2  4  *22<X2  *  1/2)* 

*  *33<V3  4  J)J  4  X44(v4  +  1)2 

*  x12(vl  4  1/2)<v2  4  l/2)  4  *13(vl  4  1/2)(v3  4  J) 


4  X1A(V1  4  ,/2)<v4  4  l)  4  *23(V2  4  1/2)(v3  4  ° 

4  *24(V2  4  ,/2)(v4  4  >>  4  *3<,(v3  4  l)l\  4  U 

4  *33  *3  4  *34  *3  2 4  4  *44  ll  4  V(J  4  ° 

♦  (*v  -  iv)K2  -  ZAJtjlj*  4  44ltK)  -  Dj  J2<J  4  l)2 

4  Djj  J(J  4  J)K2  -  »K*‘. 


(6) 
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The  rotational  constants  in  Eq.  (6)  are  further  defined  by 


Bv  -  Bg  -  Oj(Vj  +  1/2)  -  ci2(v2  +  1/2)  "  a3(v3  +  “  a4(v4  + 


Av  *  Ae  “  al(vl  +  1/2)  "  a2(v2  +  1/2)  '  a3(v3  +  °  "  a4(v4  +  1} 


Absorption  frequencies  for  Q-branch  transitions  as  deterir.ined  fror.  the 


energy 


expression  in  Eq.  (6)  and  the  selection  rules,  4v. =  1,  AJ  =  0, 


| AK (  ■  ±1,  Al3  *  AK,  are  given  by 


u  -  u  '+  2  Ta  (1  -  C,)  -  B  -  Dw  J(J  +  1)1  n  +  (a?  -  aA)n2  _  4D  n3 .  (g) 

n  vJ[_v  3vJK  J  11  K 


The  quantities  and  n  in  Eq.  (8)  are  defined  by 


W  =(jJ  -  a.  J ( J  +  1 ) 
vJ  v  3 


+  2x33(v3  +  3/2)  +  X  j  3  ( +  1/2)  +  X22(-x2  +  1  ^ 


+  X34(v4  +  1)  +  g33(213  +  1)  +  g34l4  +  A v ( 1  -  2;3)  -  Bv 


-  2Av  hh  +  'A1* J  * 


+  |k|  for  A ! K |  -  +1  R-sub  Q-branch 


- |K I  for  A  [K I  -  -1  P-sub  Q-branch 


As  disucssed  previously,  Eq.  (8)  can  be  inverted  so  that  the  quantum  label  n 

can  be  expressed  as  a  continuous  function  of  the  transition  frequency,  u. 

3 

The  cubic  term  4D^n  is  a  small  perturbation  even  for  large  values  of  n.  By 
first  ignoring  the  cubic  term,  the  resultant  quadratic  equation  for  n  can  be 
solved  to  yield 


'V>  '  «3>  -  \  '  DJKJ(J  +  ■>] 


1  J[V  '  V  ‘  Bv  -  DJK  J<J  4  »  ]  2  -  <»3  -  “!><- 


-  \J>  !/  <10’ 


/  A  ^3 

<o3  -  c3). 


The  cubic  term  can  be  introduced  as  a  perturbation  to  the  two  roots  for  n 
in  Eq.  (10).  This  is  accomplished  by  redefining  u  by 


u>  -  4D,,n” 
vj  k 


where  u',  is  substituted  for  w  ,  in  Eq.  (10).  This  process  can  be  perforate 
VJ  VJ 

iteratively  until  converged  values  for  the  roots  of  n  are  determined.  The 
root  which  is  dropped  by  approximating  Eq.  (8)  as  a  quadratic  equation 
corresponds  to  the  returning  part  of  the  P-sub  Q-branch.  Even  at  the 
highest  temperatures  of  interest  here  (-3000  K)  this  root  makes  a  negligible 
contribution  to  the  calculated  spectral  absorption  coefficients. 


The  average  spacing  of  Q-branch  transitions  about  a  given  transition 
is  given  by  dS(u >)  -  (dn/dw)-1  which  for  n  as  defined  in  Eqs.  (10)  and  (11) 
becomes 


where  the  notation,  t,  corresponds  to  the  ±  roots  of  n  in  Eq.  (10). 

The  P,  R-branch  absorption  frequencies  as  determined  by  the  energv 
expression  in  Eq .  (6)  and  the  selection  rules?6£v3  *  1 ,  £J  *  +l ,  ) 1  =  t\ 
£.1^  =  LK,  are  given  by 

*  *\*+  2[»v  -  <»3,2>  -  V2]'  <U' 

(a®  +  2Dj^.n)m^  -  4DjTt^ 
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B.  Average  Integrated  Line  Strength 

The  average  integrated  line  strength  for  Q-branch  transitions  is  giver. 
t  36-38 


S^(w)  =  S3  [(v3  +  2)/2j  exp  £  -C2  E(»/.  1,  J,  K)/tJ 


x  A  (2J  + 
JK 


D  w  £  1  -  exp  (-f2  w/T)J  j (4  i.;  T  Zv  Zr)  . 


In  Eq.  (18),  S^(w)  is  expressed  as  a  continuous  function  of  the  transition 

frequency,  w,  through  its  dependence  on  the  continuous  quantur.  nur.ber,  r.. 

defined  in  Eqs.  (10)  and  (11).  The  quantity  A  in  Eq.  (18;  corresponds  t: 

J  K 

36 

the  Honl-London  intensity  factors,  which  are 


(J  +  2  t  I K !  )  ( J  4  1  t  !  K ! ) 

(J  +  1 ) ( 2 J  +  1) 


;  LJ  *  +1  i  !k 1  =  ±] 


.(J  ±  i  t  1k!)(j  ±  Ik!) 
j(j  +  i) 


;  LJ  =  0 


!k :  =  ij 


(J  -  l  ±  Ik  I ) (j  s  IkI) 

J(2J  +  1) 


;  LJ  =  -l  i  !k  '  =  ±1 
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where  the  A„,  are  also  continuous  functions  of  u  because  of  the  expres- 
JK 

sions  for  J  and  K  in  Eqs.  (15)  and  (10).  The  vibrational  and  rotational 
partition  functions^’  ^  are  computed  from 


X  X  X  X  X  X  exp[_c2  e(  ■  °- o)/i] 

V1  v2  v3  *3  V4  *4 


Zr  -  exp  (C2  Bv/T )yj(l2Y  -i- 

V  V 


['•A 


(C2  B  /T) (1  -  \/\) 


+  !  (T/(C2  V)  ([“j  +  »JK  4  °K  ]/Av)  +  5  (T/(C2  \>)  (  ,f  .  ) 
4  2(t/(C2Bv))  (VBv)4  (T/(Vv>)  (V\) 


The  expression  for  the  integrated  line  strength  for  the  P,  R-brancr. 
p  K 

lines,  S(il),  is  almost  identical  to  the  Q-branch  line  strength  formula  in 
v 

Eq.  (18)  except  that  u/j  is  replaced  by  of  Eq.  (16). 


C.  Average  Spectral  Absorption  Coefficient  and  Line  Density 

The  average  integrated  line  strength  and  line  spacing  determine  an 

1  3  4 

average  absorption  coefficient,  which  is  expressed  as  ’  * 


D.  Effects  of  Nuclear  Spin  Statistics  on  Calculated  Band  Model  Parameters 
The  theoretical  expressions  for  the  band  model  parameters  S/d  and  1/d 
developed  in  the  preceding  sections  neglect  the  effects  that  nuclear  spin 
statistics  can  have  on  these  quantities.  An  exact  treatment  of  nuclear  spin 
would  entail  determining  the  overall  symmetry  species  of  everv  vibration- 
rotation  line  and  then  assigning  the  proper  nuclear  spin  statistical 
weight  to  each  line.  As  discussed  in  the  following  section,  there  are  far 
too  many  allowed  transitions  for  an  molecule,  like  EF^ ,  for  a'  line-by-lin 
approach  to  be  practical.  It  is  possible,  however,  to  estimate  the  average 
effects  nuclear  spin  statistics  have  on  calculated  band  model  parameters.  A 
band  average  value  of  the  absorption  coefficient  can  be  determined  fror 


<  S/d  > 


C-^ 


where  Aw  represents  an  effective  width  for  the  entire  absorption  band,  \  and 
r  represent  the  vibrational  and  rotational  quantum  states,  respectivelv,  g 
is  the  nuclear  spin  statistical  weight36  for  the  v,  r’th  transition,  c  is  the 
rotational  symmetry  number36  I  is  the  nuclear  spin36  of  the  interchanges! lc 
identical  atoms,  n  is  the  number  of  interchangeable  identical  atoms,  and  ? 
is  the  line  strength  for  the  v,  r’th  line.  It  is  important  to  realize  that 
Svr  in  Eq'  is  de*ined  neglecting  nuclear  spin  weights  entirely  (not  the 


same  as  the  case  of  I  =  0).  For  a  fundamental  band  and  its  associated  hot  bands 
it  is  well  known  both  theoretically  and  experimentally^  that  for  harmonic 
oscillators  the  total  integrated  band  strength  is  a  constant,  independent  of 
temperature,  and  it  can  be  shown  that 


v 


v 


y  s  .  s0 

i_j  vr  0 


r 


(25) 


The  implication  of  Eq .  (25)  is  that  on  average,  ^  (gvr  )/( 21  +  l)r‘J  =  1. 

Essentially,  ignoring  nuclear  spin  has  no  effect  on  S/d(_). 

A  band  average  line  density  can  be  written  as 


<  1  /d> 


(26) 


The  ratio  of  the  line  density  including  spin  effects,  <l/d>,  to  the  density 
ignoring  spin,  <l/d>’,  can  be  approximated  dv 


<  1  /d  > /  <l/d>' 


o~ 

(21  +  l)n 


i  -  1 


(27) 
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where  m  is  the  minimum  number  of  vr  states  associated  with  the  repeating 
pattern  of  nuclear  statistical  weights.  Two  simplifying  assumptions  were 
made  in  the  derivation  of  Eq.  (27)  from  Eq.  (26).  One  assumption  is  that 
for  a  fixed  vibrational  band,  v,  the  line  strength  is  a  slowly  varying  function 
of  the  rotational  indices,  r.  The  other  assumption  is  that  the  majority  of 
v,  r  bands  for  a  given  molecule  have  the  same  repeating  pattern  of  statistical 
weights.  In  general  both  these  approximations  are  reasonable  for  polyatomic 
molecules  at  room  temperature  and  above. 

The  effect  of  including  nuclear  spin  statistics  in  the  line  density  cal¬ 
culation  for  a  number  of  common  molecular  geometries  and  spin  types  is  shown, 
in  Table  II.  A  general  conclusion  that  can  be  drawn  from  Eq .  (27)  and  Table  II 
is  that  only  nuclear  spin  types  with  1=0  will  have  line  densities  substantial! 
different  from  the  line  density  computed  ignoring  the  effect  of  nuclear  spin. 

In  BF^,  ignoring  nuclear  spin  statistics  amounts  to  only  a  3”  error  in  the 
computed  line  densities. 

E.  Computational  Considerations 

A  number  of  practical  considerations  arise  in  using  a  computer  to 
evaluate  the  expressions  for  the  absorption  coefficient,  Eq.  (21),  and  the 
line  density,  Eq.  (24).  The  major  limitation  to  evaluation  of  these  band 
model  quantities  as  written  involves  the  number  of  states  that  must  be 
considered  in  the  multiple  summation  over  vibrational  states  and  rotational 
states.  For  a  MX^  molecule  six  vibrational  sums  are  necessary  (vj,  v2*  v^. 


1^.  v^,  1^)  as  well  as  additional  sums  over  the  J  and  K  rotational  quantur 
numbers  for  the  Q-branch  and  P,  R-branches  respectively.  The  number  of  states 
that  must  be  included  to  Insure  convergence  of  the  calculated  band  model 
parameters  will  scale  approximately  as  (vmax>^  J  .  For  BF^  at  high  tempera- 

g 

tures  (-3000K)  we  estimate  that  at  least  *10  quantum  states  contribute  signifi¬ 
cantly  to  the  band  model  parameters.  An  additional  multiplicative  factor  of 

2 

-10  must  also  be  included  since  the  band  model  parameters  are  evaluated  over 
a  frequency  grid,  u,  with  a  grid  spacing  comparable  to  the  spacing  between 
lines  in  a  vibration-rotation  band.  Thus,  approximately  10^  separate  evalua¬ 
tions  of  the  expressions  for  S^w)  (Eq.  (18))  and  d^(u)  (Eqs.  (12)  and  (17)) 
are  required  for  BF^  at  high  temperatures.  Although  it  is  certainly  possible 
to  perform  101®  evaluations  of  S  (u)  and  d^(u>)  on  some  of  the  larger,  faster 
computers  we  were  interested  in  a  computer  program  which  could  run  quickly  on 
the  PRIME  400  minicomputer  at  Aerodyne  Research.  A  simple  approximation  was  adopted 

which,  for  the  BF0  case,  reduced  the  number  of  evaluations  of  S  (u)  and  d  (U1 

J  v  v 

10  4 

from  -10  to  -10  .  That  is,  we  assumed  that  all  the  hot  bands  in  the  c  1 
sequence  have  the  same  normalized  band  shape.  Therefore,  once  the  absolute 
band  shape  is  known  for  a  single  band  in  the  sequence,  all  other  band  shapes 
can  be  calculated  by  scaling  to  the  appropriate  integrated  band  strength  and 
shifting  the  band  shape  to  the  appropriate  band  origin.  The  reference  band 
shape  is  calculated  from  vibrationally  average  rotational  constants.  These 
average  quantities  are  determined  by 


.■I# 1 1  rp  i 


<A  > 
v 


X  Av  (- 


C~  E(v,  0, 


v 


(26) 


where  an  analogous  expression  can  be  written  for  <B  >.  The  line  strength  of 
an  arbitrary  state,  S^,  (w),  is  related  to  the  reference  line  strength,  S  (u), 
by 

V<">  ‘  Sv("-4Vv)  [(v3+1)/(v3+1)]  («*P  [-Cj  E(v\  0,  0>/:]  / 
exp  [-C,  E(v,  0,  0) /t]  )  (7[,  -  t.v.J) 

(2  9/ 

*  (%/  [%  +  *Vv])  (  1  *  exp[‘  C2u/T])/ 

0  -  exp  C2(^  -  ^-v.v)/tJ) 
where  ^uviv  “  -  <^v<  The  line  spacings  are  related  by 


d  •  (u)  *  d  (u  -  u  ,  ) . 

V  V  V  V 


(30) 
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In  a  previous  paper^  on  linear  triatomic  band  model  parameters  several 
factors  affecting  convergence  were  discussed.  Essentially,  the  6ame  comments 
apply  here.  That  is,  the  choice  of  upper  limits  to  the  vibrational  sums  should 
be  subject  to  the  restriction  that  only  vibrational  states  whose  energy, 

E(Vj,  Vj,  v3,  1^,  v4,  1^,  0,  0),  is  less  than  an  energy  cutoff,  E  ,  should 
be  considered.  The  choice  of  Ena>,  depends  on  temperature,  where  the  rule  of 
thumb  given  previously  was 


C0  E  /T  *  10. 

2  max 


Also,  additional  terms,  constrained  to  be  positive,  were  added  to  the 
vibrational  energy  terms  in  Eq.  (6),  and  have  the  form 


Ev  *  Xin(Vj  +  1/2)3  +  x222(v2  +  1/2)3  +  X333(v3  +  l)3 


+  x444(v4  +  ’ 


(31) 


(32  i 


The  values  for  the  x,„  »  are  determined  such  that  the  first  derivative  of 

ijk  s 

the  vibrational  energy  with  quantum  number  is  always  equal  or  greater  than 

zero  for  the  maximum  quantum  numbers  considered  at  a  given  temperature  (see 

4 

Eq.  (24)  of  the  earlier  paper  for  a  prescription  to  compute  the  X^^'s^' 


V.  RESULTS  AND  DISCUSSION 


Comparisons  of  theoretical  predictions  and  experimental  determination 
of  the  absorption  coefficients  of  the  band  of  BF^  at  298,  585,  870,  1150, 

1440,  and  2400  K  are  presented  in  Figs.  6-11.  The  BF^  vibrational  and  rotation¬ 
al  molecular  constants  used  to  generate  the  theoretical  predictions  are  four.: 
in  Table  III.  Most  of  the  BF^  molecular  constants  vere  determined  experimental!' 
from  high  resolution  room  temperature  spectral  measurements.  Those  molecular 

constants  not  obtainable  from  room  temperature  data  have  been  estimated  by 
20 

Ginn  et.  al.,  who  considered  the  effect  of  anharmonic  potential  energy  terms 
in  the  BF^  intramolecular  potential  energy  function  in  the  prediction  of  these 
constants.  For  the  vibrational  mode  the  molecular  constants  which  art  nr: 
determined  from  low  temperature  spectral  measurements  are  >^3,  >33.  >3  . 

g3^i  and  834-  At  room  temperature  852  of  the  BF3  molecules  are  in  the  ground 
vibrational  state,  and  most  of  the  remaining  157  are  in  the  v  =  1  level.  T':.t 
maximum  displacement  of  the  hot  bands  associated  with  the  doubly  degenerate 
*=  1  level  is  at  most  -7  cm  *  (X34  +  834)  which,  given  the  much  smaller  inter.si 
of  this  band  compared  to  the  fundamental,  means  that  the  perturbatiom  this 
has  on  the  bandshape  is  relatively  small  and  therefore  difficult  to  measure 
accurately.  In  order  to  more  accurately  determine  these  molecular  constants 
it  is  necessary  to  perform  high  temperature  spectral  measurements  where  the 
hot  band  transitions  are  sufficiently  intense  relative  to  the  rotational  envelope 
of  the  fundamental  band.  The  sensitivity  of  the  predicted  band  contour  at 
2400  K  to  changes  in  the  molecular  constants  >^3  and  X34  is  shown  in  Fig.  11. 
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TABLE  III.  HBF3  Molecular  Constants20,  22 *  36  used  In  Band 
Model  Calculations. 


Molecular  Constant^ 

Value  (cm-*) 

Molecular  Constant 

Value  (err.  ^ ) 

“1 

940. 

A 

e 

0.17635 

“2 

699. 

B 

e 

0.35:' 

w(b) 

u-3 

1486.5 

A 

4  x  1C'" 

j 

*"4 

481. 

Q1 

A 

-3  x  10~4 

q2 

Ml 

-0.99 

A 

7.5  x  10'4 

x22 

1.31 

A 

Q4 

1  x  1C'~4 

i  x33 

-12.28 

B 
a , 

8  x  10'" 

1 

' 

1  x44 

0.39 

B 

-1.5  x  1C" 

X12 

-3.23 

B 

U3 

—  X 

1.5  x  1C 

X1 3 

-4.10 

B 
a  , 

-4.7  x  1C" ■' 

4 

X14 

-8.12 

0.795 

v3 

x23 

-11.0 

-0.811 

s4 

X24 

0.92 

D, 

4.1  x  1C"' 

J 

X34 

2.0 

n 

-6.7  x  10*' 

JK 

833 

9.66 

n 

2.7  x  10*' 

K 

g34 

4.67 

844 

-0.20 

(a) 

(b) 


Observed  fundamentals  corrected  for  anharmonicity  (x^)  to  give 
harmonic  frequencies  Uj,  <*>3,  and  J 

For  *°BF3  ■  1537.5  cm  *.  The  ^B/^°B  isotope  ratio  used  in  the 

calculations  was  4.10. 
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Initial  estimates  of  X23  “  “2.11  and  X34  e  3.8  cm  1  were  used  (see  curve  "A"  of 

20 

Fig.  11)  which  were  the  values  predicted  by  Ginn  et.  al.,  who  included  the 
effect  of  anharmonicity  in  the  BF^  intramolecular  potential  energy  surface  or. 
the  prediction  of  vibrational  and  rotational  molecular  constants.  By  adjust¬ 
ing  these  constants  to  x23  *  and  X34  *  2  cm  a  good  fit  to  the  flame 

spectra  at  2400  K  was  attained.  Since  it  is  clear  that  the  high  temperature 
data  can  be  fit  well  with  only  two  parameters,  v„,  and  Xo,»  it  is  not  possible 

to  simultaneously  determine  x^.  X33,  833.  and  g^ .  The  assumed  values  of 

20 

*13’  *33’  ^33’  ^34*  were  predicted  by  Ginn  et.  al .  One  expects  that 

X23  and  X34  will  have  the  greatest  influence  on  the  high  temperature  bandshape 

because  and  are  the  lowest  energy  modes  in  BF^-  At  temperatures  below 

1500  K  (flow  tube  data)  the  predicted  contours  were  much  less  sensitive  tc 

variations  in  x23  and  X34;  although  the  set  of  constants  determined  by  fittir.; 

the  2400  K  data  did  produce  a  significantly  better  fit  than  the  initial  estimates . 

There  are  two  areas  of  disagreement  between  the  data  and  the  theoretical 

predictions  for  the  temperatures  below  2400K.  The  first  of  these  is  the  fact 

that  the  data  in  the  center  of  the  band  is  generally  higher  than  the  predictions 

and  the  experimental  widths  of  the  band  are  slightly  narrower  than  the  prediction 

at  585,  870,  1150,  and  1440K.  In  addition,  for  all  of  the  flow  reactor  data  it 

was  not  possible  to  obtain  a  consistent  value  for  the  oand  strength.  The  band 

-2  -1 

strengths  derived  from  the  data  ranged  from  310C  to  5900  cm.  amagat  .  At  a 
given  temperature,  the  band  strength  measured  on  different  days  showed  a  simi¬ 
lar  variability.  However,  the  measured  band  shapes,  normalized  to  a  common 
band  strength,  were  repeatable.  Both  of  these  disagreements  apparently  arise 
from  the  same  problem  -  the  recirculation  of  infrared  active  species  in  the 


White  cell  optics  tubes.  The  observed  effects  could  be  caused  by  relatively 
little  recirculation  since  the  optical  path  through  the  recirculating  gases  is 
much  longer  than  that  through  the  main  flow  and  due  to  the  relatively  long  resi¬ 
dence  time  of  species  in  a  recirculating  flow. 

It  can  be  seen  from  Table  1  that  the  flow  reactor  was  always  operated 
at  Reynolds  numbers  which  should  be  high  enough,  as  indicated  by  the  flow 
simulation  experiment,  to  preclude  recirculation  of  infrared  active  species 
in  the  optics  tubes.  In  addition,  a  check  experiment  was  conducted  at 
870K  in  which  the  Reynolds  number  was  doubled,  with  no  resulting  significant 
change  in  the  measured  absorption.  This  seems  to  indicate  that  the  recir¬ 
culation  is  not  due  to  operation  of  the  flow  reactor  at  a  "marginal"  Reynolds 

number.  However,  in  the  flow  simulation  experiment  the  shield  flow  was 

43 

observed  to  recirculate.  Thus,  it  is  likely  that  the  recirculation  may  bt 
due  to  diffusion  of  infrared  active  species  into  the  shield  flow.  This 
phenomenon  was  not  observed  in  the  flow  simulation  experiment  because  the 
residence  time  of  a  volume  of  fluid  in  the  analysis  region  scales  linearly 
with  the  Reynolds  number,  but  the  diffusion  length  scales  as  the  square  root 
of  the  residence  time.  Therefore,  for  a  given  Reynolds  number  the  diffusion 
length  In  the  gas  will  be  greater  than  the  diffusion  length  in  the  flow 
simulation  experiment. 

For  future  experiments  the  nozzle  on  the  end  of  the  main  flow  tube  will 
have  its  inside  diameter  reduced  from  4.9  cm  to  4.2  cm.  This  will  reduce  the 
ratio  of  main  nozzle  area  to  (concentric)  shield  nozzle  area  from  2:1  to 
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approximately  1:1  (with  a  concurrent  increase  in  the  shield  flow).  This 
increase  in  the  width  of  the  shield  flow  in  the  analysis  region  should  eliminate 
the  recirculation  of  infrared  active  species  from  the  main  flow  in  the  optics 
tubes. 

The  integrated  band  strength  of  the  band  of  BF^  has  been  previously 

23  —2  —I 

determined  by  McKean  at  room  temperature  to  be  3290  ±  300  cm  amagat 

This  can  be  compared  to  the  2400K  value  determined  by  the  flame  measurements 

of  2800  ±  230  cm  amagat  .  It  has  been  shown  by  Breeze  et.  al.,  that  in  the 

harmonic  oscillator  approximation  the  integrated  band  strength  is  independent 

52 

of  temperature.  However,  Crawford,  and  Dinsmore  have  demonstrated  that  vher. 

the  effects  of  mechanical  and  electrical  anharmonicity  are  considered,  the 

bandstrength  is  no  longer  temperature  independent.  The  band  model  formulatic:. 

used  in  the  theoretical  predictions  of  the  BF^  band  contours  partially  includes 

the  effect  of  mechanical  anharmonicity  by  using  the  vibrational  anharmonic 

molecular  constants  in  the  evaluation  of  the  vibrational  partition  functier. 

-2-1 

and  excited  state  populations.  Using  an  input  value  of  3290  cm  amagat 

for  the  room  temperature  bandstrength  the  model  predicts  bandstrengths  of 

-2  -1  -2  -1 
3170  cm  amagat  at  1500  K  and  2940  cm  amagat  at  2500  K.  This  suggests 

that  the  difference  in  the  experimental  bandstrength  determinations  may  be 

real  and  can  be  explained  by  the  effects  of  mechanical  and  electrical  anhar.oni- 

city  on  the  temperature  dependence  of  integrated  bandstrengths. 

A  complete  set  of  theoretical  predictions  for  the  S/d  and  1/d  band  model 

parameters  at  300,  600,  900,  1500,  2000,  2500,  and  3000  K  can  be  found  in 
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Figs.  12-17.  The  extremely  large  line  densities  above  600  K  preclude  the 
experimental  determination  of  the  1/d  parameter  because  y(l/d)  >>  1,  where  y  is 
the  linevidth  of  a  single  transition,  even  when  y  is  Doppler  limited.  At  room 
temperature  and  below,  1/d  is  measureable  at  low  pressure  (Doppler  limit).  An 
attempt  to  measure  1/d  at  room  temperature  was  made;  however,  the  results  were 
inconclusive  due  to  the  recirculation  problem  discussed  above. 


3-41 


ii  niiir  nlttfcfiisrrriiHiii  <n  ■’iiil'  furiliiiil 


ACKNOWLEDGEMENTS 


We  wish  to  thank  R.  Brown  for  his  major  role  in  the  construction  of  the 
flow  reactor  and  in  his  assistance  in  obtaining  the  data  and  J.  Moore  for  his 
dedicated  work  in  obtaining  the  burner  spectra.  We  also  wish  to  acknowledge 
valuable  discussions  with  C.  Kolb,  M.  Camac,  F.  Kaufman,  H.  Wolfhard, 

S.  Zakanycz,  F.  Simmons,  G.  Lindquist,  and  C.  Arnold. 

The  work  at  AR1  was  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  Contract  F49620-77-C-0075.  The  work  at  ER1M  was  sponsored 
by  che  Defense  Advanced  Research  Projects  Agency  and  monitored  by  MI COM . 


3-42 


REFERENCES 


1.  C.B.  Ludwig,  W.  Malkmus,  J.E.  Reardon,  and  J.A.L.  Thomson,  "Handbook 
of  Infrared  Radiation  from  Combustion  Gases, "NASA  SP-3080,  (Scientific 
and  Technical  Information  Office,  Washington,  D.C.,  1973). 

2.  S.S.  Penner,  R.C.  Sepucha,  and  J.E.  Lowder,  J.  Quant.  Spectrosc.  Radiat. 
Trans.  10,  1001  (1970). 

3.  W.  Malkmus,  J.  Opt.  Soc.  Am.  53,  951  (1963). 

4.  L.S.  Bernstein,  J.  Quant.  Spectrosc.  Radiat.  Trans.  23,  157  (1980). 

5.  L.S.  Bernstein,  D.C.  Robertson,  and  J.A.  Conant,  J.  Quant.  Spectrosc. 
Radiat.  Trans.  23,  169  (1980). 

6.  W.M.  Elsasser,  Phys.  Rev.  54,  126  (1938). 

7.  G.N.  Plass,  J.  Opt.  Soc.  Am.  48,  690  (1958). 

8.  C.C.  Ferisso,  C.B.  Ludwig,  and  L.  Acton,  J.  Opt.  Soc.  Am.  56,  17]  (1966). 

9.  Y.  Ben-Aryeh,  Appl.  Opt.  6,  1049  (1967). 

10.  W.0.  Davies,  J.  Opt.  Soc.  Am.  54,  467  (1964). 

11.  U.P.  Openheim  and  Y.  Ben-Aryeh,  J.  Opt.  Soc.  Am.  53,  344  (1963). 

12.  T.  Kunitomo  and  M.  Osumi,  J.  Quant.  Spectrosc.  Radiat.  Trans.  15,  345  (397 

13.  L.D.  Gray,  Appl.  Opt.  4,  1494  (1965). 

14.  S.G.W.  Ginn,  J.K.  Kenney,  and  J.  Overend,  J.  Chem.  Phys.  48,  1571  (1968). 

15.  S.G.W.  Ginn,  C.W.  Brown,  J.K.  Kenney,  and  J.  Overend,  J.  Mol.  Spectrv.  28, 

509  (1968). 

16.  C.W.  Brown  and  J.  Overend,  Can.  J.  Phys.  ^6,  977  (1968). 

17..  C.W.  Brown  and  J.  Overend,  Spectrochim.  Acta  25A,  1535  (1969). 

18.  S.N.  Dreska  and  K.N.  Rao,  J.  Mol.  Spectrv.  _1_8»  404  (1965). 

19.  S.G.W.  Ginn,  D.  Johansen,  and  J.  Overend,  J.  Mol.  Spectry.  36>  448  (1970). 

20.  S.G.W.  Ginn,  S.  Reichman,  and  J.  Overend,  Spectrochim.  Acta  26A  291  (1970) 

21.  J.L.  Duncan,  J.  Mol.  Spectry.  J_3,  338  (1964). 

3-43 


AD-A089  743  AERODYNE  RESEARCH  INC  BEDFORD  MA  CENTER  FOR  CHEMICAL— ETC  F/9  7/4 

QUANTITATIVE  INFRARED  SPECTROSCOPIC  MEASUREMENTS  OF  BORON  FLUOR— ETC  (U> 
AU6  80  ME  6ERSH*  C  E  KOLB#  M  S  ZAHNISER  F49620-77-$-0075 

UNCLASSIFIED  ARl-RR-225  _  AF0SR-TR-80-0942  NL 


- ■  wtm*  iimiuwp  iihkb4I,i 


22.  J.L.  Duncan,  J.  Mol,  Spectry,  22 ,  247  (1967). 

23.  D.C.  McKean,  J.  Chem.  Phys.  24,  1002  (1956), 

24.  I.W.  Levin  and  S.  Abramowitz,  J.  Chem,  Phys,  £3,  4213  (1965), 

25.  D.T.  Hawkins,  L.S.  Bernstein,  W.E.  Falconer,  and  W.  Klemperer,  "Binary 

Fluorides:  Free  Molecular  Structures  and  Force  Fields,  A  Bibliography 
(1957  -  1975)."  (Plenum  Pub.  Corp.,  New  York,  N.Y.,  1976). 

26.  C.R.  Bailey,  J.B.  Hale,  and  J.K.  Thompson,  J.  Chem.  Phys.  5,  274  (1937). 

27.  C.R.  Bailey,  J.B.  Hale,  and  J.W.  Thompson,  Proc.  Roy.  Soc.  A  161,  107  (1937). 

28.  P.A.  Freedman  and  W.J.  Jones,  J.  Mol.  Spectrosc.  54^,  182  (1975). 

29.  D.M.  Gage  and  E.F.  Barker,  J.  Chem.  Phys.  7,  455  (1939). 

30.  S.  Konaka,  Y.  Murata,  K.  Kuchitsu,  and  Y.  Morino,  Bull.  Chem.  Soc.  Jap. 

39,  1134  (1966). 

31.  K.  Kuchitsu  and  S.  Kinaka,  J.  Chem.  Phys.  45,  4342  (1966). 

32.  F.N.  Masri,  J.  Mol.  Spectrosc.  43,  168  (1972). 

33.  G.M.  Sheldrick,  J.  Mol.  Spectrosc.  20,  295  (1966). 

34.  J.  Vanderryn,  J.  Chem.  Phys.  30,  331  (1959). 

35.  D.M.  Yost,  D.  Devault,  T.F.  Anderson,  and  E.N.  Lassettre,  J,  Chem.  Phvs. 

6,  424  (1938). 

36.  G.  Herzberg,  "Infrared  and  Raman  Spectra,"  (Van  Nostrand  Reinhold,  New  York, 
1945). 

37.  W.F.  Edgell  and  R.E.  Moynihan,  J.  Chem.  Phys.  45,  1205  (1966). 

38.  S.L.  Gerhard  and  D.M.  Dennison,  Phys.  Rev.  43,  197  (1933). 

39.  R.M.  Badger  and  L.R.  Zumwalt,  J.  Chem.  Phys.  b,  711  (1938). 

40.  D.M.  Dennison,  Phys.  Rev.  318  (1926). 

41.  G.  Placzek  and  E.  Teller,  Z.  Physik  til.,  209  (1933). 


3-44 


42.  T.H.  Edwards,  J.  Opt.  Soc.  Am.  51^  98  (1961). 

43.  M.E.  Gersh  and  C.E.  Kolb,  In  "Characterization  of  High  Temperature 
Vapors  and  Gases"  (Edited  by  J.W.  Hastie),  NBS  Special  Publication 
561  (1979).  p.  679. 

44.  A.  Fontijn,  S.C.  Kurzius,  J.J.  Houghton,  and  J.A.  Emerson,  Rev.  Sci. 

Instr.  43,  726  (1972). 

45.  E.E.  Ferguson,  F.C.  Fehsenfeld,  and  A.  Schmeltekopf ,  Adv.  At.  Mol.  Phvs. 

5,  1  (1969). 

46.  R.A.  Oetjen,  C.  Kao,  and  H.M.  Randall,  Rev.  Sci.  Instr.  1_3,  515  (1942). 

47.  F.S.  Simmons,  "Spectroscopic  Pyrometry  of  Gases,  Flames  and  Plasmas," 

ISA  Transactions,  2,  (1963),  p.  168. 

48.  A.G.  Gaydon  and  H.G.  Wolfhard,  "Flames,  Their  Structure,  Radiation  and 
Temperature,"  (Chapman  and  Hall,  Ltd.,  London,  1969),  2nd  Edition. 

49.  S.  Gordon  and  B.J.  McBride,  Report  No.  NASA-SP-273  (National  Aeronautics 
and  Space  Administration,  Washington,  D.C.,  1971). 

50.  The  terms  involving  the  centrifugal  distortion  constants  were  derived 
using  the  formulation  in  N.  Davidson,  "Statistical  Mechanics"  (McGraw- 
Hill  Co.,  New  York,  1962),  p.  172. 

51.  J.C.  Breeze,  C.C.  Ferriso,  C.B.  Ludwig,  and  W.  Malkmus,  J.  Quant.  Spectrosc. 
Radiat.  Trans.  402  (1965). 

52.  B.L.  Crawford  and  H.L.  Dinsmore,  J.  Chem.  Phys.  2®,  983  (1950). 


3-45 


Figure  Captions 

1.  Cross-sectional  view  of  the  Aerodyne  Research,  Inc,,  heated  flow  reactor, 

Shield  gas  flow  (A),  main  gas  flow  (B),  7.3  cm  alumina  flow  tube  (C) , 

water  cooled  vacuum  jacket  (D),  refractory  insulation  (E),  six-arm  cross 
to  support  flow  tubes  (F),  entry  and  exit  points  of  infrared  radiation 
(G),  gas  flow  to  pumps  (H),  mass  spectrometer  sampling  orifice  (I), 
exhaust  port  to  pump  (J),  white  cell  mirror  (1  of  2)  (K),  optics  tube 

purging  gas  flow  (1  of  4)  (L),  Kanthal  heating  elements  (M) ,  and 

multihole  mullite  heater  array  (N) . 

2.  Cross-sectional  view  of  flow  reactor  analysis  region.  Main  flow  (A), 
purge  flow  (B),  mass  spectrometer  sampling  orifice  (C),  six-arm  support 
cross  (D),  and  refractory  insulation  (E). 

3.  External  optics  of  flow  reactor.  Flat  mirrors  (M2,  M3,  M4,  M7,  M9), 
spherical  mirrors  (Ml,  M8),  off-axis  paraboloid  mirrors  (M5,  M6), 

HgCdTe  detector,  and  Spex  Minimate  monochromator. 

4.  Optical  schematic  of  Enviromental  Research  Insitute  of  Michigan  double 
beam  emission-absorption  spectrometer.  Shutters  (B1  -  B4),  choppers 
(Cl,  C2),  blackbody  <BB),  spectrometer  (SP),  detector  (D) ,  CaF  windows 
(W),  and  grating  (G). 

5.  Equilibrium  calculations  of  boron  species  combustion  products  for  BF^ 

flame  measurements  at  2400  K.  ' 

6.  Comparison  of  theoretical  and  experimental  (flow  tube)  determinations  of 
the  BF^  band  absorption  coefficients  at  298K. 

7.  Comparison  of  theoretical  and  experimental  (flow  tube)  determinations  of 
the  BF^  band  absorption  coefficients  at  585K. 

8.  Comparison  of  theoretical  and  experimental  (flow  tube)  determinations  of 
the  BF^  band  absorption  coefficients  at  870K. 

9.  Comparison  of  theoretical  and  experimental  (flow  tube)  determinations  of 
the  BF^  v-  band  absorption  coefficients  at  1150K. 

10.  Comparison  of  theoretical  and  experimental  (flow  tube)  determinations  of 
the  BF^  band  absorption  coefficients  at  1440K. 

11.  Comparison  of  theoretical  and  experimental  (flame)  determinations  of 
the  BF^  band  absorption  coefficients  at  2400K.  Initial  estimates20 
of  *  -2.11  and  x-ia  *  3*8  cm”1  (A)  and  adjusted  values  of  x->,  *  -11. 
and  X34  -  2.  cm"!  (B}7 
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Theoretical  prediction  of  the  BF,  v,  band  absorption  coefficients  at 
300  and  600K. 

Theoretical  prediction  of  the  BF.  v,  band  absorption  coefficients  at 
900,  1200,  and  1500K.  J  J 

Theoretical  prediction  of  the  BF.  v.  band  absorption  coefficients  at 
2000,  2500,  and  3000K.  J 

Theoretical  prediction  of  the  BF,  v.  band  line  density  parameter,  1/d, 
at  300  and  600K. 

Theoretical  prediction  of  the  BF.  v.  band  line  density  parameter,  1/d, 
at  900,  1200,  and  1500K. 

Theoretical  prediction  of  the  BF,  v.  band  line  density  parameter,  1/d, 
at  2000,  2500,  and  3000K.  *  J 
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INPUT  BF3  MOLE  FRACTION 

Figure  5 
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it.  MEASUREMENT  OF  THE  VIBRATIONAL  BAND  STRENGTH  OF  nBF 
USING  A  TUNABLE  DIODE  LASER 


by 

Mark  S.  Zahniser  and  Michael  E.  Gersh 

ABSTRACT 

The  infrared  band  strength  for  ^BF  is  measured  by  line  absorption  from 
the  P(ll)  transition  at  1343.5  cm  *  using  a  tunable  diode  laser  and  multi¬ 
pass  absorption  cell  coupled  to  a  discharge-flow  apparatus.  The  BF  number 
density  is  simultaneously  determined  by  ultraviolet  absorption  in  the 

unresolved  band  head  of  the  A  II  X  transition  at  196  nm.  The  result- 

-2  -1 

ing  value  for  the  band  strength  is  650  +  200  cm  (STP  atm)  which  is  in 
good  agreement  with  existing  theoretical  calculations. 
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I.  INTRODUCTION 


The  electronic  spectrum  of  boron  monofluoride  has  been  extensively 

12  2 

studied  in  the  ultraviolet  in  both  emission  '  and  absorption  and  its  pure 

3 

rotational  spectrum  has  been  observed  by  microwave  spectroscopy.  There 
have  been,  however,  no  observations  in  the  infrared  of  its  vibrational- 
rotational  spectrum.  This  is  often  the  case  for  transient  species  due  to 
the  difficulties  of  producing  them  in  sufficient  quantity  for  detection  by 
conventional  infrared  techniques. 

There  are  even  fewer  IR  vibrational  band  strength  measurements,  which 
require  an  independent  determination  of  absolute  concentration,  for  short¬ 
lived  species.  Such  measurements,  however,  have  been  recently  obtained  for 
the  CS-0  radical”4,  ^  and  the  N0+  ion^  using  tunable  diode  lasers. 

In  this  work  we  report  the  direct  observation  of  BF  using  a  diode  laser 
at  7.3  urn  to  measure  the  integrated  absorption  of  a  single  vibrational-rotational 
line.  The  line  strength  is  determined  by  measuring  the  BF  concentration  simul¬ 
taneously  by  UV  absorption  using  the  band  head  of  the  A  *11  X  *1  (0,0)  transi¬ 
tion  at  196  nm.  The  results  of  this  study  provide  a  method  of  monitoring  BF 
concentrations  as  required  for  laboratory  kinetic  studies.  They  are  also  of 
interest  for  comparison  with  ab  iniffo  calculations  of  molecular  properties 
of  simple  diatomic  molecules,  particularly  the  dipole  moment  derivative  with 
respect  to  internuclear  distances. ^  Finally,  the  measured  band  strength  and 
the  known  spectroscopic  constants  may  be  used  to  produce  an  accurate  IR  radia- 

g 

tion  band  model. 
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II.  EXPERIMENTAL 


Measurement  of  the  IR  absorption  coefficient  for  transient  species 
requires:  (1)  a  stable  source  of  the  molecule,  (2)  a  sensitive  detection 
system,  and  (3)  a  well  defined  optical  path  which  must  be  long  enough  (several 
meters)  to  provide  sufficient  optical  depth  at  low  concentrations  and  yet  must 
be  traversed  by  the  active  species  in  a  time  short  enough  (several  milliseconds) 
so  that  the  concentration  remains  constant  over  the  path.  Ideally,  the  path 
should  not  be  bounded  by  surfaces  such  as  windows  or  mirrors  on  which  the 
active  species  may  recombine  and  thereby  create  a  concentration  gradient  along 
the  optical  path. 

These  requirements  are  met  in  the  apparatus  shown  in  Figs.  1  and  2  which 

9 

combines  a  discharge  flow  system  with  a  multipass  absorption  "White"  cell 
aligned  so  that  the  optical  path  is  perpendicular  to  the  direction  of  flow. 

The  alumina  flow  tube  is  7.3  cm  id  and  1  m  long,  but  is  reduced  in  diameter 
to  4  cm  where  it  joins  the  detection  region.  At  this  point  a  concentric  flow 
of  argon  is  added  to  the  main  flow  and  adjusted  so  that  the  velocities  of  the 
two  flows  are  equal.  This  results  in  an  aerodynamic  confinement  of  the  flow 
tube  gases  in  the  central  portion  of  the  detection  region  and  has  the  advantages 
of:  (1)  maximizing  the  column  density  of  IR  active  species  within  the  absorption 
path,  (2)  minimizing  their  traversal  time  across  the  folded  path,  and  (3)  pre¬ 
venting  labile  species  from  coming  into  contact  with  the  White  cell  mirrors  or 
other  surfaces  in  the  detection  region.  A  more  detailed  description  of  the 
apparatus  and  of  the  conditions  for  which  the  flow  is  effectively  confined  is 
given  elsewhere.*® 
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BF  is  produced  in  a  microwave  discharge  (2450  MHz,  100  W)  of  BF^  (Matheson, 

99.5%  pure)  and  helium  in  a  1.0  cm  id  alumina  tube  which  is  connected  to  the 

main  flow  tube  by  a  1.4  cm  id,  60  cm  long,  alumina  tube.  Maximum  BF  concentra- 
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tions  are  obtained  with  relatively  large  flow  rates  of  helium,  150  STP  cm  s  , 

3  -1 

and  much  smaller  flow  rates  of  BF^,  1  STP  cm  s  .  BF^  is  observed  by  IR  absorp¬ 
tion  at  7.3  pm  and  found  to  be  less  than  1%  dissociated  by  the  discharge.  Argon, 

3  -1 

100  STP  cm  s  ,  is  added  to  the  main  flow  tube  to  give  a  total  pressure  of  2 

torr  and  a  linear  velocity  of  2500  at  s  \  BF  concentrations  as  high  as 
12  -3 

2  x  10  molecules  cm  are  obtainable  under  these  conditions.  All  experi¬ 
ments  in  this  study  are  done  at  300K. 

A.  UV  Absorption  Measurements 

BF  concentrations  are  determined  by  UV  absorption  in  the  unresolved  band 
head  of  the  A  X  *1  transition  at  196  nm.  The  light  source  is  a  high  pres¬ 

sure  xenon  arc  lamp  (Varian  Model  VIX-500UV)  which  has  a  total  continuum  output 
of  -2  mW  nm  *  from  195  to  200  nm.  The  light  is  collimated  by  an  internal  para¬ 
bolic  reflector,  and  baffles  on  both  ends  of  the  absorption  cell  define  a  0.5  cm 
diameter  beam  which  passes  through  the  center  of  the  flow  field  and  the  IR  absorp¬ 
tion  axis  and  is  perpendicular  to  both.  The  beam  is  focused  by  a  quartz  lens 
onto  the  entrance  slit  of  a  1  m  monochromator  (McPherson  Model  2051).  The 
grating  (1200  lines/mm,  300  nm  blaze)  is  driven  by  a  digital  stepping  motor  and 
is  used  in  second  order  with  10  pm  slits  which  give  a  resolution  of  0.007  nm. 

A  200  nm  interference  filter  with  20  nm  bandpass  (Acton  Research)  transmits  less 
than  0.01%  of  the  first  order  radiation  at  400  nm.  The  phototube  (Hamamatsu  R763P) 
is  operated  in  analog  mode.  Signal  noise  due  to  lamp  intensity  fluctuations  is 
reduced  by  splitting  the  beam  before  it  enters  the  absorption  cell,  monitoring 


it  with  a  second  photomultiplier  through  a  200  nm  interference  filter,  and 
comparing  the  two  signals  using  a  differential  amplifier. 

The  fractional  UV  absorption  is  measured  by  turning  the  microwave  discharge 

on  and  off  while  the  monochromator  wavelength  is  at  a  fixed  position  within  the 

l  1 

R-branch  band  head  of  the  A  II  •*-  X  £  transition.  The  data  analysis  requires  the 
relative  monochromator  transmission  at  each  spectral  line  position  within  the  band¬ 
pass,  which  depends  on  the  slit  function,  G^,  and  the  absolute  wavelength  setting 
of  the  monochromator.  The  slit  function  is  determined  by  scanning  through  the 
Hg  lines  from  a  low  pressure  pen  lamp  (Oriel,  Model  6037)  at  19A.1  nm  in  second 
order  and  312.56  nm  in  first  order  and  by  calibrating  the  scan  rate  using  the 
latter  line  and  the  adjacent  313.15  nm  doublet.  The  spectral  area  under  the  slit 
is  determined  by  integrating  over  the  transmitted  line  shape  with  a  polar  plani- 
meter  which  gives  /  G\  dX  *  0.0091  +  0.0002  nm  in  second  order.  The  slit  function 
is  empirically  approximated  in  a  functional  form  by  a  Gaussian  core  with  Lorent- 
zian  tails  for  the  purpose  of  calculating  the  relative  transmission  for  each  BF 
spectral  line.  Line  positions  for  1JBF  are  calculated  from  the  molecular  constants 
given  by  Caton  and  Douglas. ^  Line  positions  for  the  10BF  isotope  are  shifted 
0.7  cm  *  to  the  red  and  are  calculated  according  to  Herzberg.^  The  spectral 
slit  encompasses  14  lines,  R(8)  through  R(21),  for  each  isotope.  The  band  head 
is  shaded  towards  the  red  and  lies  between  the  atmospheric  absorption  lines  of 
the  3,0  Schumann  Runge  band  of  Oj.  Although  the  monochromator  is  purged  with 
N2,  there  is  still  sufficient  O2  in  the  optical  path  to  provide  a  reference 
for  determining  the  absolute  monochromator  setting  by  scanning  from  the  R(19) 

I  _  1  to 

and  P ( 1 7)  O2  lines  at  51129.49  cm  and  51121.74  cm  ,  respectively,  through 
the  BF  band  head  at  51111.26  cm  The  spectrum  is  recorded  with  a  strip  chart 
recorder,  and  the  scan  is  stopped  at  51110.0  cm  *  which  is  slightly  to  the  red 
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of  the  absorption  maximum.  The  uncertainty  in  the  resulting  monochromator 
setting  is  ±  0.4  cm  1  which  corresponds  to  an  uncertainty  of  ±5%  in  deter¬ 
mining  from  the  convolution  of  the  slit  function  with  the  absorption 

lines. 

The  fractional  absorption  from  a  continuum  source  which  has  a  constant 
intensity  over  the  monochromator  bandpass  is 


where  N  is  the  number  density  of  the  isotope,  I,  in  a  rotational  state,  J; 

I  ,  J 

£  is  the  pathlength;  0Q(J)  *  2Sj/hv^  V  n2  is  the  absorption  cross  section 

at  line  center,  vq,  for  the  Doppler-broadened  rotational  line;  u  =  2  V  £n2 

(v  -  v0)/£''>d  i-s  the  reduced  frequency  in  terms  of  the  Doppler  width,  Av^  =  0.116 

cm  S  and  is  the  monochromator  slit  transmission  function.  The  integrated 

absorption  cross  section,  S^,  is  calculated  from  the  radiatiive  lifetime,  t, 
i 

of  the  A  II  electronic  state, 


1 

J  0  2  g"  t  ^oo  2J  +  1 

8tt  c  v  6 


(2) 


where  vq  is  the  transition  frequency  in  cm  *,  g'  *  2  and  g"  *  1  are  the  electro¬ 
nic  degeneracies  of  *11  and  *1  states,  respectively,  q  is  the  Franck-Condon 
factor  for  the  0,0  vlbronic  band,  and ^  ■  (J  +  2)/2  is  the  Honl-London  factor*1 
for  the  R-branch  transitions  from  the  ground  state  which  has  been  normalized 
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so  that  I  =  2J  +  1,  where  J  and  J*  refer  to  the  lower  and  upper  states, 

J’  J 

1  7 

respectively.  The  lifetime,  t  =  2.8  i  0.3  ns  has  been  measured  by  Hesser 

using  the  electron  excitation-phase  shift  technique.  Franck-Condon  factors 

14  15 

have  been  calculated  by  Wentink  and  Spindler  and  by  Mishra  and  Khanna  who 
obtain  qQo  *  0.79  and  0.78,  respectively.  The  dependence  of  the  electronic 
transition  probability  on  the  variation  in  the  r-centroid  is  small  enough*"’ 
to  be  neglected  in  determining  S^.  The  fractional  population  in  a  rotational 
level  for  a  given  isotope  is 


¥■  -  s  Bo 2<J  -  •** 


[- 


B  J(J  +  1)  hc/kT 
o 


] 


(3) 
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where  the  rotational  constant  is  Bq  e  1.507  cm  .  The  fractional  isotopic 

abundances  for  ^BF  and  10BF  are  0.80^  and  0.198,  respectively . *^ 

The  total  BF  number  density  is  determined  from  a  plot  of  fractional 

absorption  versus  £bfJ  obtained  by  numerically  integrating  Eq.  (1).  The 

maximum  concentration  used  in  these  experiments  is  £  BF J  =  1.7  x  10*^  molecules 
-3 

cm  which  corresponds  to  a  fractional  absorption  of  7.3%  and  a  maximum  optical 

depth  at  any  frequency  of  INjjOoJl~0.3  with  the  path  length  £  =  4.3  cm.  The 
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detection  limit  of  2  x  10  molecules  cm  is  determined  by  instabilities  in 
the  arc  lamp  which  limits  the  minimum  detectable  change  in  AI/I  to  -1%  with 


a  signal/noise  ratio  of  unity  and  a  10  s  integration  time  constant. 


B.  IR  Absorption  Measurements 

The  infrared  detection  system  uses  a  Pb-salt  diode  laser  (Laser  Analytics) 
contained  in  a  closed  cycle  helium  refrigeration  unit  with  temperature  regula¬ 
tion.  A  KRS-5  lens  focuses  the  radiation  onto  the  entrance  slit  of  a  0.25  m 
monochromator  (Spex  Minimate).  A  150  line/mm  grating  with  250  pm  slits  provides 
sufficient  resolution,  1  cm  \  to  separate  laser  frequency  modes  and  gives  a 

coarse  wavelength  calibration.  Precise  wavelength  calibration  is  obtained  by 

17  18 

comparison  with  known  methane  absorption  lines  ’  from  a  gas  reference  cell 
placed  in  the  beam.  The  beam  is  mechanically  chopped  at  800  Hz  for  direct  absorp¬ 
tion  experiments,  and  brought  to  a  focus  at  the  plane  of  the  front  mirror  of  the 
White  cell.  The  White  cell  mirrors  have  a  43  cm  radius  of  curvature  and  are 
adjusted  for  40  passes  as  a  compromise  between  increased  path  length  and  decreased 
beam  intensity  due  to  reflectance  losses.  The  number  of  passes  is  determined  by 
counting  the  number  of  spots  on  the  front  mirror,  which  are  formed  by  a  helium- 
neon  laser  which  is  co-aligned  with  the  IR  beam. 

The  IR  beam  is  focused  onto  a  liquid  nitrogen  cooled,  HgCdTe  detector 
(Santa  Barbara  Industries)  whose  output  signal  goes  to  a  phase  sensitive  ampli¬ 
fier  referenced  to  the  mechanical  chopper.  Aosorption  spectra  are  displayed 
using  an  XY  recorder  where  the  X-axis  is  driven  by  the  diode  current  which  is 
proportional  to  the  laser  frequency.  The  frequency  tuning  rate  of  9.1  t  0.3 
cm  */A  is  determined  by  scanning  through  the  CH^  reference  lines  and  is  independ¬ 
ently  confirmed  by  the  frige  spacing  from  an  accidental  86  cm  etalon  formed  by 
a  slight  misalignment  of  the  back  mirrors  of  the  White  cell  which  results  in 
multiple  reflections  of  one  edge  of  the  expanded  beam. 
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The  particular  laser  mode  used  in  these  experiments  has  a  continuous 
tuning  range  from  1341.5  cm  *  to  1344.1  cm  * .  Modes  in  the  spectral  regions 
on  either  side  of  this  range  are  continuous  over  a  much  smaller  range  before 
mode  hopping  occurs,  and  therefore,  the  P(ll)  line  of  the  fundamental  band 
of  ^BF  at  1343.5  cm  ^  was  chosen  for  this  study.  The  laser  is  operated  at 
a  temperature  of  60 K  and  a  current  of  1.3  A;  conditions  which  often  result 
in  multimodes  too  closely  spaced  to  be  separated  by  the  monochromator.  The 
mode  used  in  this  study,  however,  was  shown  to  be  singular  by  adding  enough 
CH^  to  the  absorption  cell  to  completely  attenuate  the  laser  intensity  in  the 
region  of  the  BF  line. 

Absorption  measurements  are  made  by  alternately  scanning  the  laser  through 
the  absorption  line  with  the  microwave  discharge  on  and  then  off  as  showTi  in 
Fig.  3.  The  spectra  are  then  overlaid  and  the  difference  is  traced  with  a 
planimeter  to  obtain  the  integrated  absorption  intensity.  This  procedure  is 
preferred  to  taking  absorption  measurements  only  at  line  center  since  it  mini¬ 
mizes  the  uncertainty  due  to  accidental  etalon  fringes  which  are  superimposed 
on  the  baseline,  and  it  eliminates  the  need  to  assume  lineshape  functions  for 
both  the  absorber  and  the  laser  source  at  these  loptical  depths.  The  lat¬ 
ter  is  an  important  consideration  since  rather  wide  variability  has  been  reported, 

-3  -1 

with  laser  line  widths  as  large  as  1.4  x  10  cm  ,  when  using  closed  cycle  cryo¬ 
stats  rather  than  liquid  helium  dewars  for  cooling.  The  measured  temperature 
tuning  rate  of  the  diode  in  the  present  study  is  0.8  cm  *  K  *  which,  combined 

with  the  manufacturer's  specification  for  cryostat  temperature  stability  of 

-3  -3  -1 

x  1  x  10  K  ,  gives  an  expected  line  width  on  the  order  of  10  cm  compared  to 


* 

s 
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the  Doppler  width  for  BF  of  3  x  10  cm  .  Measurements  are  taken  under 

optically  thin  conditions  so  the  maximum  fractional  absorption  in  these 

experiments  is  less  than  4%  at  line  center.  The  laser  source  is  sufficiently 

-3 

stable  to  detect  fractional  absorptions  as  low  as  10  directly  with  a  Is  time 

constant,  which  corresponds  to  a  detection  limit  for  BF  of  -4  x  lO1^  molecules 

-3 

cm  .  Although  greater  sensitivity  is  possible  by  modulating  the  source  and 

20 

using  second  derivative  detection,  the  direct  absorption  method  is  better 
suited  for  line  strength  determinations.  The  problem  of  a  sloping  baseline 
which  is  usually  characteristic  when  frequency  tuning  diode  lasers  is  minimized 
here  by  adjusting  the  monochromator  setting  so  that  the  change  in  intensity 
due  to  current  tuning  the  laser  is  offset  by  the  change  in  the  transmission 
function  of  the  monochromator. 

C.  Optical  Path  Length  Determinations 

Since  the  UV  and  IR  absorptions  are  measured  simultaneously  on  the  same 
volume  element,  only  the  path  length  ratio  is  needed  to  interpret  the 

results.  This  ratio  is  just  the  integral  number  of  passes  of  the  IR  beam  in  the 
White  cell  provided  that  both  UV  and  IR  beams  are  small  compared  to  the  flow  dim¬ 
ensions  and  pass  through  the  center  of  a  uniform  flowfield.  To  confirm  these 
conditions  the  ratio  is  independently  verified  by  making  absorption  measurements 
on  stable  gases,  NO  in  the  UV  and  CH^  in  the  IR,  using  the  identical  flow  con¬ 
ditions  and  geometry  as  in  the  BF  experiments.  The  effective  path  length  is 
obtained  by  comparing  the  observed  absorption  for  a  measured  mixing  ratio  of 
absorber  in  the  carrier  gas  under  dynamic  conditions  to  the  absorption  observed 
with  static  conditions  when  the  flow  tube  is  isolated  from  the  pumping  system 
and  filled  with  a  measured  partial  pressure  of  the  stable  absorber.  The  path 
lengths  in  the  static  case  are  the  directly  measureable  distances  between  the 


outer  windows  in  the  UV,  56  cm,  and  between  the  White  cell  mirrors, 


40  x  43  cm,  in  the  IR. 

ip  -l17’18 

The  R(6)E1  line  in  the  fundamental  of  12CH^  at  1342.94  cm 

is  used  for  the  IR  measurements.  The  in  1°.  for  intensities  at  line  center 

1 

is  plotted  versus  [CH0  for  both  the  static  and  dynamic  conditions  and  for 

fractional  absorptions  from  10  to  60%.  The  effective  IR  path  length  for  the 

dynamic  case  is  calculated  as  the  ratio  of  the  slopes  of  these  plots  times  the 

static  path  length  and  is  £.  *  173  ±  5  cm,  or  4.3  cm  for  a  single  pass.  Although 

the  path  length  determination  using  this  ratio  method  does  not  depend  on  the 

absolute  value  of  the  CH^  line  absorption  coefficient,  the  corresponding  value 

for  the  integrated  line  strength  from  this  experiment,  assuming  a  Doppler  line 

-20  -1  -2 

shape,  is  S  -  2.8  x  10  cm  /molecule  •  cm  .  This  is  in  good  agreement  with 

-20  21 
the  value  of  S  «  2.7  x  10  derived  using  the  line  parameters  of  Chedin  et  al . 

22 

and  line  absorption  measurements  of  Restelli  at  al.  on  other  lines,  and  it 

23  24 

is  within  the  range  of  band  strength  measurements  by  several  other  authors.  ’ 

The  UV  path  length  measurement  uses  line  absorption  by  the  Rj2(10.5)  and 

(^OO. 5)  degenerate  lines  of  the  0,0  vibronic  band  of  the  A  2Z  *■  X2r  transition 
-1  25 

in  NO  at  44,135  cm  which  is  fully  resolved  by  the  1  m  monochromator  in 

second  order  with  10  pm  slits.  The  path  length  is  calculated  from  the  fractional 
absorption  according  to  the  relationship 


41 
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F(N.a  H) 
J  o 

/Gv  dv 


(4) 


where  Nj/N  is  the  fractional  population  in  the  lower  rotational  state  and  the 
00  _ 

term  F(N.  0  l)  •  Tn  (-N.  0  l)  / (n+1) ! V n+i  accounts  for  the  finite  optical 
j  o  n“0  J  o 

26 

depth  at  line  center.  The  maximum  optical  depth  in  these  experiments  is  1.3 
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which  corresponds  to  a  fractional  absorption  of  7%.  The  ratio  of  absorption 

measurements  under  static  and  dynamic  conditions  gives  a  value  of  = 

3.8  ±  0.8  cm  where  the  larger  uncertainty  compared  to  the  IR  measurement  is 

due  to  the  greater  uncertainty  in  the  fractional  absorption  measurement.  The 

UV  path  length  may  also  be  calculated  directly  from  Eq.  (4)  using  the  integrated 

—  18  2  —1  -1 

line  absorption  coefficient  of  (Nj/N)Sj  ■  4.0  x  10  cm  molecule  cm  obtained 

27 

from  the  0,0  band  oscillator  strength.  This  results  in  a  value  of  = 

4.4  ±  0.8  cm  from  absorption  measurements  taken  with  the  dynamic  conditions  only. 

The  combined  UV  and  IR  static  versus  dynamic  ratio  measurements  give  a  value 
^IR^UV  '  ^  *  9  and  confirm  the  integral  path  length  ratio  of  40,  which  is  used 
in  the  BF  data  analysis,  within  experimental  uncertainty.  These  measurements 
also  confirm  the  effectiveness  of  the  aerodynamic  confinement  of  the  flow  tube 
gases  to  the  center  of  the  detection  region  and  imply  that  there  is  negligible 
recirculation  of  the  absorbing  species  in  the  side  arms.  This  is  further  con¬ 
firmed  by  visual  observation  in  the  detection  region  of  an  afterglow  produced 
by  passing  air  through  the  microwave  discharge  with  the  same  flow  conditions 
that  were  used  for  the  BF  absorption  measurements.  The  afterglow  was  seen  to 
be  well  confined,  with  no  apparent  recirculation. 
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III.  RESULTS  AND  DISCUSSION 


The  BF  line  frequency,  1343.47  +  0.02  cm  1 ,  is  determined  by  interpolation 

J 

between  two  sets  of  methane  lines  as  shown  in  Fig.  3.  The  relative  CH^  line 

18 

positions  are  in  good  agreement  with  those  of  Blatherwick  et  al.  which  have 

an  absolute  accuracy  of  ±0.005  cm  *  .  The  line  position  calculated1*  for  the 

J1BF  P(ll)  line  is  1343.38  cm  1  using  the  molecular  constants  obtained  from  the 
12  3 

electronic  spectra  ’  and  microwave  spectra  which  are  tabulated  by  Huber  and 
28 

Herzberg.  The  rotational  constants  obtained  from  the  UV  spectra  by  Ca ton 
2 

and  Douglas  are  in  excellent  agreement  with  those  obtained  from  the  microwave 
3  -1 

spectra,  and  the  0.1  cm  discrepancy  between  the  calculated  and  observed  line 
position  is  most  probably  associated  with  the  uncertainty  in  the  vibrational 
constraints,  and  o^x^  which  are  obtained  from  the  band  origins  in  the 
UV  spectra.*  The  observed  line  position  corresponds  to  a  value  of  vq  *  1378.73 
for  1,0  vibrational  band. 

In  addition  to  the  close  agreement  between  expected  and  observed  line 
position,  further  evidence  that  the  line  is  unambiguously  due  to  1JBF  is  pro¬ 
vided  by  the  observation  that  no  other  absorption  lines  due  to  the 
microwave  discharge  are  apparent  in  the  2.5  cm  1  region  covered  by  the  laser 
mode.  The  spacing  between  adjacent  J1BF  lines  is  2Bq  *  3.0  cm  1  and  neither 
the  P(10)  and  P ( 1 2)  lines  could  be  located  with  this  diode  due  to  frequency 
gaps  in  the  tuning  range.  The  10BF  P(21)  line  is  separated  from  the  nBF  P ( 1 1 ) 
by  only  0.3  cm  1 ,  but  it  is  an  order  of  magnitude  less  intense  due  to  the  higher 
J  value  and  is  not  observed.  Line  positions  for  the  2,1  vibrational  band  are 


also  in  this  region  and  would  be  observed  if  there  were  a  significant  unrelaxed 
non-Boltzmann  distribution  of  vibrational  levels  formed  in  the  microwave  dis¬ 
charge.  The  other  possible  species,  BF2,  may  be  eliminated  on  line  strength 
and  line  spacing  arguments  since  the  larger  partition  function  for  the  triatomic 
species  greatly  decreases  the  intensity  for  a  single  rotational  transition  com¬ 
pared  to  a  diatomic  species  with  the  same  band  strength,  and  more  than  one  line 
would  be  apparent  in  a  2.5  cm  '  interval  due  to  the  larger  moment  of  inertia. 

There  is  also  thermodynamic  evidence  against  the  existence  of  large  quanti- 

29 

ties  of  BFj  in  this  system. 

The  line  width  in  Fig.  3  is  30%  greater  than  the  Doppler  width  at  300 K 
of  3.0  x  10  cm  .  Several  other  effects  may  contribute  to  the  linewidth 
including:  (1)  pressure  (Lorentz)  broadening,  (2)  unresolved  hyperfine  split¬ 

ting,  and  (3)  instrumental  broadening  due  to  the  laser  linewidth.  Pressure 
broadening  coefficients  are  typically  on  the  order  of  0.05  -  0.1  cm  *  atm 
(half-width  at  half  maximum),  for  N2  and  02  as  third  bodies  and  are  expected 
to  be  less  with  Ar  as  the  collision  partner  as  in  this  case.  The  experimental 
lineshape  may  be  fitted  with  a  Voigt  profile,  but  this  requires  a  broadening 
coefficient  of  0.2  cm-1  atm-1  which  is  unrealistically  large  for  Ar.  Although 

it  is  possible  that  the  broadening  coefficient  for  BF^  collisions  could  be 

15  -3 

significantly  greater,  the  maximum  BF^  concentration  is  <10  molecules  cm  which 
is  too  low  to  account  for  this  effect.  The  magnitude  of  the  hyperfine  interaction 
contribution  to  the  linewidth  may  be  estimated  from  the  value  for  eq  Q  '  -4.5  MHz 


for  nB19F  from  microwave  spectra-5  and  is  negligible  compared  to  the  Doppler 

width.  The  most  probable  cause  of  the  larger  than  Doppler  linewidth  is, 

therefore,  instrumental  broadening.  Although  we  have  no  direct  measurement 

-3  -1 

of  the  laser  linewidth,  a  value  on  the  order  of  1  x  10  cm  is  not  unreason¬ 
able,  as  mentioned  previously,  and  could  account  for  the  difference.  This  is 
further  supported  by  observations  of  the  CH^  line  profiles  which  also  are  also 
larger  than  Doppler.  It  is,  or  course,  quite  possible  that  both  Lorentz 
broadening  and  instrumental  broadening  contribute  to  the  observed  lineshape. 
However,  for  low  optical  depths  and  integrated  absorption  measurements  it  is 
not  necessary  to  know  the  exact  line  profile  to  determine  the  line  strength. 

The  line  strength  is  calculated  from  the  absorption  data  given  in 
Table  1  using  the  relationship 


where  the  integrated  absorption  intensity  is  calculated  from  the  fractional 
absorption  at  line  center  times  the  ratio  of  area.  A,  to  height,  h,  of  the 
absorption  profile,  /(AI/Iq)  dv  *  (AIMA^/Io)  A/h;  the  product  | £ BfJ 
is  the  total  BF  column  density  from  the  UV  absorption  measurements;  the 
factor  0.802  is  the  fractional  isotopic  abundance  of  J1BF;  and  *uy/£jR  “  1/AO 
is  the  path  length  ratio.  The  above  is  exact  only  in  the  limit  of  low  optical 
depth.  The  approximations  introduced  by  not  Integrating  over  the  exact  line- 
shapes  for  both  laser  Intensity  and  absorption  coefficient  contribute  no  more 


3 


than  2%  to  the  overall  uncertainty  in  these  experiments,  where  the  maximum 

optical  depth  is  0.05  assuming  a  doppler  lineshape  and  a  laser  linewidth  of 
-3  -1 

10  cm  .  The  average  value  for  five  absorption  measurements  at  each  of  two 

-19  2  -1 

different  BF  concentrations  is  =  (7.2  +  2)  x  10  cm  molecule 

cm-*,  at  300 K,  which  corresponds  to  a  fundamental  band  strength  of  SgANI)  = 

—2  —1 

650  1  200  cm  (STP  atm)  .  The  error  limits  represent  the  overall  uncertainty, 
2o,  where  the  major  contributions  are:  the  random  errors  associated  with  the 
UV  absorption  measurements,  ±  20%  ;  the  combined  accuracy  of  the  UV  lifetime 
measurement  and  the  Franck-Condon  factors  used  to  determine  the  UV  absorption 
coefficient,  ±15%  ;  the  measured,  integrated  IR  line  profile,  ±10%  ;  and  the 
slit  function  convolution  for  the  UV  band  head  absorption,  ±5%. 

The  band  strength  of  a  diatomic  molecule  is  directly  related  to  the  dipole 
moment  derivative  with  respect  to  internuclear  distance  by 


S 


-vl  ur 

3c\  M 


(6) 


where  N  is  Loschmidt’s  number  and  m  is  the  reduced  mass.  For  dp/dr  in 
o  r 

units  of  Debye/X,  and  the  reduced  mass  in  amu,  this  gives  S  =  188  j'dp/drJ'Vm^ 

-2  -1 

cm  (STP  atm)  .  The  results  of  several  theoretical  calculations  of  dy/dr  are 
compared  with  the  experimental  value  of  4.9  +  0.8  Debye/X  in  Table  II.  The 
recent  configuration  Interaction  calculations  (Cl)  of  Kurtz  and  Jordan  repre¬ 
sent  the  most  sophisticated  treatment  and  are  in  good  agreement  with  the  experi¬ 
mental  value.  The  agreement  is  also  within  the  experimental  error  for  both  the 

30  31 

SCF  calculations  of  Nesbet  and  of  Huo  .  The  semiempirical  CNDO  calculations 
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which  were  performed  with  a  program  suppled  by  the  Quantum  Chemistry  Program 
32 

Exchange  give  a  result  which  is  larger  than  either  the  experimental  value 

or  the  ab  initio  calculation;  however,  the  agreement  is  reasonably  good, 

33  34 

as  is  typical  of  CNDO  calculations  of  dipole  moment  derivatives.  ' 

The  large  band  strength  for  BF  makes  it  feasible  to  use  infared  absorp¬ 
tion  as  a  diagnostic  in  reaction  kinetic  studies.  Fractional  absorptions  as 
small  as  10  ^  are  readily  detectable  with  tunable  diode  laser  sources  using 

9 

derivative  techniques  which  corresponds  to  a  detection  limit  of  <10  molecules 

_3 

cm  using  the  experimental  configuration  in  this  work.  This  sensitivity  is 

35 

comparable  to  UV  resonance  fluorescence  which  has  been  previously  used  in 
3(> 

BF  kinetic  studies,  and  has  the  additional  advantage  compared  to  the  resonance 
fluorescence  technique  of  providing  absolute  concentration  measurements  without 
further  calibration. 


i 

!■ 
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TABLE  I.  BF  Absorption  Data 


TABLE  II.  Comparison  of  Theoretical  Calculations  of  Dipole 
Moment  Derivative  for  15BF  with  Present  Experimental  Results 
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Figure  Captions 

1.  Top  view  of  flow  tube  with  detail  of  analysis  region  showing 
confinement  of  optically  active  gases  by  the  shield  flow.  Flow 
tube  length  is  not  to  scale.  White  cell  (A),  microwave  discharge 
cavity  (B) . 

2.  Cross  sectional  view  of  analysis  region  with  infrared  absorption 
and  ultraviolet  absorption  diagnostics.  Elements  external  to 
vacuum  jacket  not  shown  to  scale.  Vacuum  jacket  (A),  six  arm 
cross  piece  (B) ,  white  cell  mirrors  (C) »  windows  (D) ,  arc  lamp  (E) , 
beam  splitter  (F),  interference  filters  (G),  photomultipliers  (H) , 
mirrors  (I) ,  quartz  lens  (J) ,  choppers  (K) ,  monochromators  (L) , 

KRS5  lens  (M) ,  diode  laser  (N) ,  HgCdTe  detector  (0).  IR  path 
within  White  cell  shows  maximum  divergence  of  laser  beam.  The 
central  dashed  circle  indicates  the  boundary  of  the  optically 
active  flow  tube  gases. 

3.  Diode  laser  scan  through  P(ll)  Line  of  ^BF  with  reference 
lines.  Insert  is  expanded  10  x  in  intensity  and  shows  scans 
with  and  without  the  microwave  discharge.  The  dashed  line 

is  with  discharge  off  sumperimposed  on  the  trace  with  dis¬ 
charge  on.  BF  concentration  is  1.0  x  10*2  molecules  cm-3 
and  path  length  is  172  cm.  The  CH4  reference  line  positions 
from  Ref.  18  are  shown  as  arrows  above  the  frequency  scale. 

CH4  is  added  to  the  flow  tube  at  the  beginning  of  the  scan, 
turned  off  at  label  A,  and  on  again  at  label  B.  The  lower 
frequency  lines  are  100  x  more  intense  than  the  higher  fre¬ 
quency  triplet  and  are  completely  absorbed  at  line  center 
which  indicates  the  singularity  of  the  laser  mode.  The 
fringes  in  the  baseline  are  due  to  the  White  cell  mirrors 
which  produce  an  etalon  effect  with  free  spectral  range  of 
0.0058  cm"1. 
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